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I .  INTRODUCTION 


This  report  presents  the  significant  results  of  a  two  and 
one  half  year  research  program  conducted  "  y  A.R.A.P,  under 
contract  number  DAAD05-76-C-0757 .  Tt_  primary  objective  of 
the  program  was  to  gain  an  understanding  of  the  underlying 
physics  of  the  impact  process  and  to  apply  this  knowledge  in 
the  development  of  a  simple  model  to  describe  the  penetrator- 
target  interaction. 

The  program  evolved  in  three  distinct  phases.  Phase  1 
dealt  with  the  characterization  of  target  materials.  Phase  2 
dealt  with  the  development  of  a  long  rod  penetrator  model. 

The  final  phase  concerned  the  measurement  of  the  partitioning 
of  energy  during  fracture  of  a  penetrator  specimen. 

Interim  reports  were  Issued  at  the  end  of  Phase  1  and 
Phase  2  (Refs.  1  and  2).  Some  of  the  work  described  in  these 
reports  was  incomplete,  and  some  of  the  results  have  had  to 
be  revised  in  the  light  of  subsequent  developments.  It  is 
the  purpose  of  this  report  to  tie  together  the  loose  ends  from 
the  first  two  phases  and  to  present  for  the  first  time  the 
results  of  Phase  3  of  the  research  program. 

Prior  to  the  inception  of  this  program,  a  useful  tool  for 
the  study  of  the  impact  process  was  under  development  at 
A.R.A.P.  This  tool,  called  the  Integral  Theory  of  Impact, 
sought  to  describe  the  Impact  process  in  as  simple  terms  as 
possible  without  sacrificing  the  essential  physics  which  were 
thought  to  describe  the  problem.  Rather  than  worry  about  the 
microscopic  details  of  the  interaction  of  two  materials 
during  Impact,  it  was  decided  to  model  the  Interaction  in  a 
global  or  integral  sense.  In  this  way,  a  set  of  ordinary 
rather  than  partial  differential  equations  was  obtained  - 
an  attractive  simplification  from  both  cost  and  time  view¬ 
points.  Since  the  equations  used  in  the  theory  retain  the 
essential  physics  of  the  impact  process  in  an  integral  sense, 


1.  Donaldson,  Coleman  duP.,  Contiliano,  Ross  M. ,  and  McDonough, 
Thomas  B. :  The  Qualification  of  Target  Materials  Using  the 
Integral  Theory  of  Impact.  A.R.A.P.  Rept.  No.  295, 
Aeronautical  Research  Associates  of  Princeton,  Inc. , 

December  1976. 

2.  Swanson,  Claude  V.  and  Donaldson,  Coleman  duP. :  Applica¬ 
tion  of  the  Integral.  Impact  Theory  to  Modeling  Long-Rod 
Penetrators.  A.R.A.P.  Rept.  No.  333,  Aeronautical  Research 
Associates  of  Princeton,  Inc.,  March  1978. 


it  was  felt  that  the  important  features  of  such  processes 
could  be  exhibited  in  spite  of  the  simplicity  of  the  model. 

Thus,  the  integral  theory  was  born  in  an  attempt  to  bridge 
the  gap  between  the  complex  and  costly  multi-element  codes 
and  the  purely  empirical  models. 

It  was  not  intended  that  the  integral  theory  replace 
either  the  large  codes  or  the  experimentalist.  Because  of 
its  simplicity,  the  model  introduces  a  degree  of  economy 
which  makes  it  reasonable  to  conduct  parametric  studies  and 
observe  trends  rather  than  single  points.  In  this  way,  the 
integral  theory  can  be  used  to  guide  experimental  programs, 
interpret  results,  and  Berve  as  a  screening  tool  for  cases 
which  reauire  the  details  which  are  only  available  in  the 
large  codes. 

Early  studies  (Ref.  3)  of  the  impact  process  showed  that 
for  the  purpose  of  computing  target  response  during  impact 
it  was  necessary  to  determine  at  least  two  characteristic 
quantities,  in  addition  to  density,  for  any  target  material. 

One  of  these,  denoted  by  E*p  and  called  the  hydrodynamic 

mode  energy,  represents  the  amount  of  energy  required  to  put 
a  well  defined  mass  of  target  material  in  a  hydrodynamic 
state.  The  other  quantity,  denoted  by  E*g  ,  represents 

the  elastic  energy  absorbed  by  the  same  mass  of  target  mater¬ 
ial  during  impact.  Phase  1  of  this  program  consisted  of  an 
experimental  and  theoretical  evaluation  of  E*p  and  E*e 

for  a  wide  spectrum  of  target  materials.  Reference  1  summa¬ 
rized  the  experimental  program  and  evaluation  of  E*p  and 

E,Vfi  and  the  development  of  equations  to  predict  these  two 

quantities  from  fundamental  material  properties.  More  recently, 
additional  materials  have  been  qualified  and  coefficients  in 
the  theoretical  equations  have  been  modified  as  more  data 
have  become  available.  In  addition,  the  correlations  which 
were  reported  in  Ref.  1  have  been  revised  to  include  a 
better  estimate  for  the  drag  coefficient  which  appears  in 
the  equations.  The  results  of  all  qualification  tests  and 
the  present  best  estimates  for  the  theoretical  value  of  FVf  ,  i.e., 
the  sum  of  E,v  and  E*  ,  are  described  in  Chapters  II  and  III. 


3.  Donaldson,  Coleman  duP. ,  McDonough,  Thomas  B.,  and 

Contiliano,  Ross  M. :  Application  of  the  Integral  Impact 
Theory  to  the  Design  of  Specialized  Tactical  Ordnance. 
A.R.A.P.  Rept.  No.  279,  Aeronautical  Research  Associates 
of  Princeton,  Inc.,  May  1976. 
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During  the  second  phase  of  this  program,  attention  was 
directed  toward  the  design  of  rod  penetrators.  The  concepts 
of  the  integral  theory  were  applied  to  the  development  of  a 
long  rod  penetrator  model.  A  two-element  model  was  developed 
which  incorporates  the  physics  of  the  deforming  single  cell 
model  developed  several  years  ago  by  A.R.A.P.  (Ref.  4)  and 
couples  this  cell  to  a  rigid  shaft.  This  model  was  described 
in  Ref.  2.  However,  the  code  which  was  developed  contained 
several  bugs  and  was  difficult  to  utilize.  In  order  to 
eliminate  these  problems,  the  equations  of  the  model  were 
recast,  and  the  code  was  rewritten.  The  revised  model  and 
code  are  described  in  this  report. 

The  rod  model  requires  the  specification  of  a  property, 

*  which  represents  the  non-recoverable  or  dissipated 

energy  required  to  hydrodynamicize  a  unit  mass  of  the  rod 
material  as  the  tip  of  the  rod  is  forced  to  flow  out  of  the 
path  of  the  relatively  undamaged  shaft  of  the  penetrator. 

In  essence,  is  to  the  rod  what  E*p  is  to  the  target. 

The  property  E*^  can  be  obtained  by  correlating  experimental 

data  with  calculations  using  the  two-cell  model.  Alternatively, 
an  experimental  technique  was  devised  by  A.R.A.P.  to  measure 
E,^  directly.  The  results  of  these  experiments  are  also 

described  in  this  report. 

In  what  follows,  a  summary  of  the  entire  target  materials 
qualification  program  is  provided  in  Chapter  II.  For  complete¬ 
ness,  the  test  data  are  provided  in  Appendix  A,  and  the  deduced 
values  of  E*  are  shown  in  Appendix  B.  The  theoretical  model 
for  E.(f  based  on  fundamental  material  properties  is  provided 
in  Chapter  III  together  with  the  present  estimates  for  the 
theoretical  E*  of  many  materials.  Chapter  IV  contains 
a  description  of  the  long  rod  penetrator  model.  A  user’s 
guide  for  the  numerical  code  is  provided  in  Appendix  C.  Chap¬ 
ter  V  describes  the  E^^  experiments,  and  Chapter  VI  contains 

the  conclusions  drawn  from  this  program. 


4.  Donaldson,  Coleman  duP.,  Contiliano,  Ross  M. ,  and  McDonough, 
Thomas  B. :  A  Study  of  Water  Drop  Displacement  and  Deforma¬ 
tion  in  Aerodynamic  Shock  Layers.  A.R.A.P.  Rept.  No.  265, 
Aeronautical  Research  Associates  of  Princeton,  Inc., 

March  1978. 
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II.  TARGET  MATERIAL  QUALIFICATION 

In  this  chapter,  the  method  by  which  the  impact  proper¬ 
ties  E,t  and  E,tg  of  a  target  material  are  obtained  will 

be  described.  Briefly,  the  procedure  consists  of  conducting 
a  series  of  impact  tests  using  nondeforming  projectiles  and 
relatively  thick  target  samples.  By  eliminating  the  deforma¬ 
tion  of  the  projectile  and  backface  effects  in  the  target, 
the  equations  of  motion  are  simplified  and  an  impact  can  be 
completely  specified  given  the  density  and  geometry  of  the 
projectile  and  target  plus  the  material  property  of  the 

target.  In  the  next  section,  the  Integral  Theory  of  Impact 
will  be  developed  for  a  nondeforming  sphere  impacting  a 
semi-infinite  target  to  show  how  E*p  and  E*e  are  evaluated. 

In  the  subsequent  section,  the  results  of  an  experimental 
program  to  characterize  a  wide  range  of  target  materials  will 
be  summarized. 

A.  Rigid  Sphere  Impacting  Semi-Infinite  Target 

1.  Decelerating  Force  -  Consider  the  penetration  of  a 
semi-infinite  target  by  a  nondeforming  sphere  of  radius  R 
and  mass  nip  which  is  traveling  at  a  velocity  VCq  normal 

to  the  target.  At  any  instant  during  the  penetration,  the 
projectile,  which  has  traveled  a  distance  z  ,  has  a  velocity 
Vc  and  is  being  decelerated  by  a  force  F  .  The  nomenclature 
is  depicted  schematically  in  Figure  1.  Conservation  of 
momentum  requires  that 


ar(Vc)-  -F  ■ 


(1) 


If  both  sides  of  Eq.  (1)  are  multiplied  by  V£  ,  the 

result  is  the  equation  for  the  rate  of  change  of  the  center 
of  mass  kinetic  energy 


d 

dt 


(2) 


The  product  FVC  is  the  rate  at  which  the  projectile  loses 

kinetic  energy  or  the  rate  at  which  work  is  done  on  the  target 
by  the  projectile.  If  W.  is  the  total  energy  of  the  target 
material,  then 


4 


d 

ai 


w. 


fv. 


(3) 


or 
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(A) 


which  is  Just  a  restatement  of  conservation  of  total  energy. 

The  energy  absorbed  by  the  target  is  distributed  between 
kinetic  energy  of  the  target  material  Kt  and  a  quantity  Ut 

which  represents  nonkinetic  forms  of  energy  such  as  plastic 
work  and  elastic  energy 


d  Wt  d  Kc  d  U 
dt  dt  dt 


(5) 


As  the  projectile  penetrates,  the  target  volume  swept  out  per 
unit  time  at  angle  4>  on  the  sphere  is  given  by 


dv 

3E 


2TtR4sin  $  dd  V. 


cos 


VcA 


(6) 


o 

where  2 ttR1  Bin  d  cos  <f>  d<f>  represents  the  elemental  area  normal 
to  the  velocity  vector  Vc  at  an  angle  <t>  from  the  nose  and 

<(>  is  the  submergence  angle, 
s 

The  rate  of  change  of  kinetic  energy  can  be  written  as 


d  Kt  dv 
dt  Pt  dt 


PtVcA 


CDVc 


(7) 


and  the  change  of  nonkinetic  energy  can  be  written  as 


<8) 


d  tL  dv 

— -  -  »t  -  <E*> 
dt  u  dt 


pt  Vc  A  E*  > 


where  pt  is  the  target  mass  density  and  Cp  represents  the 

instantaneous  drag  coefficient  of  the  projectile.  It  will 
be  shown  that  CD  is  related  to  the  Newtonian  drag  coefficient 
of  the  sphere.  u 

The  term  E*  in  Eq.  (8)  consists  of  two  parts:  a  constant 
part  which  represents  plastic  work  and  a  term  E*e  which 

represents  elastic  energy  and  is  significant  for  small  penetra¬ 
tions.  The  product  pt  E*p  ,  which  accounts  for  plastic  flow 

work  and  has  the  units  of  pressure  is  analogous  to  the  Brinell 
hardness  of  the  target  material.  However,  it  models  the  flow 
strength  of  the  material  at  the  strain  rate  of  impact  rather 
than  at  the  quasi-static  strain  rate  of  a  Brinell  test.  In 
reality,  it  is  an  adiabatic  rather  than  an  isothermal  measure¬ 
ment  of  the  vield  strength  because  the  time  scale  for  shear 
heating  in  the  material  is  much  faster  than  the  time  scale 
for  thermal  relaxation.  The  quantity  pt  E*p  is  therefore 

called  the  "adiabatic  hardness"  of  the  material.  It  will  be 
shown  later  that  pt  E*p  is  roughly  constant  as  a  function 

of  velocity  for  each  material.  In  Ref.  1,  a  formula  was 
derived  for  E*p  based  on  Brinell  hardness,  melting  tempera¬ 
ture,  and  specific  heat.  A  summary  of  this  work  is  given  in 
Chapter  III. 

The  quantity  E*e  represents  the  pressure  produced  by 

elastic  compression  of  the  target  material.  A  formula  for 
E*e  was  also  derived  in  Ref.  1,  based  on  the  elastic  modulus, 

Brinell  hardness,  and  density  of  the  target  material.  At  the 
end  of  this  section,  the  method  by  which  E*  is  evaluated 
from  impact  tests  will  be  described. 

Comparison  of  Eqs.  (5),  (7),  and  (8)  with  Eq.  (3)  3hows 
i hat  the  decelerating  force  can  be  written  as 


F 


Pt  A 


(s 


(9) 
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Based  on  Newtonian  flow  considerations  but  allowing  some 
freedom  for  comparison  with  experimental  data,  it  is  possible 
to  write  the  force  in  the  following  form 


2  r 8  / 

’  0  2  \ 

F  -  2tt p c  R^  1  sin  <j>  cos  <t>  f 

cosVE*p  +  E*ej 

d*  .  (10) 


For  this  analysis,  the  quantity  a  which  is  related  to  the 
Newtonian  drag  coefficient  is  constant  and  has  the  value 
a  *  1  .  Equation  (10)  can  be  integrated  to  yield 


(11) 


where  $  is  given  by 
8 


4>  ■  cos_1(l  -  z/R)  z  <  R 

o  " 


«  v/2 


z  >  R 


(12) 


Note  that  the  integration  of  Eq.  (10)  can  only  be  done  in  this 
manner  if  E*e  is  independent  of  $  .  Thus,  the  term  EVfe 

which  appears  in  Eq.  (11)  is  an  average  value  of  the  elastic 

energy  for  the  submergence  angle  .  It  is  not  the  local 

s 

value  of  the  elastic  energy  per  unit  mass  absorbed  by  the 
target.  More  on  this  point  later. 

With  (11)  for  the  force,  the  momentum  equation  (1)  becomes 


d 

dt 


3  1  P, 


4Rp, 


(13) 


In  Eq.  (13),  the  definition  for  proj ectile  mass ,  mp  -  (4/3)  ("R^pp)  , 
has  been  used.  The  penetration  is  governed  by 


8 


dz 

3t 


(14) 


Equations  (13)  and  (14)  together  with  Eq .  (12)  can  be  inte¬ 
grated  simultaneously  to  yield  z(t)  and  Vc(t)  .  The 

integration  proceeds  from  t  -  0  to  t  -  tf  ,  the  time  at 

which  Vc  is  reduced  to  0  . 

2,  The  Newtonian  drag  coefficient  -  It  was  noted  above 
that  the  drag  force  is  based  on  the  Newtonian  flow  approxima¬ 
tion.  In  this  approximation,  the  target  material  is  assumed 
to  flow  in  a  very  thin  layer  along  the  surface  of  the  pro¬ 
jectile.  As  a  small  volume  of  target  material  enters  this 
layer,  it  suffers  an  abrupt  change  in  direction  as  it  comes 
in  contact  with  the  projectile.  In  the  Newtonian  model,  the 
normal  component  of  velocity  is  destroyed  while  the  tangential 
component  is  conserved  as  the  control  volume  enters  the  shear 
layer.  It  is  the  destruction  of  the  normal  component  of 
velocity  which  Increases  the  pressure  at  the  surface  of  the 
projectile.  The  momentum  transferred  by  each  volume  of 
material  is 


d  dv 

— (momentum)  »  p*.  —  V„  cos  <j>  .  (15) 

dt  c  dt  c 


The  component  of  this  momentum  in  the  axial  direction  is  the 
drag  force.  The  contribution  to  the  total  drag  force  at  the 
angle  $  is  given  by 


Fdrag 


ov  2 

P,.  —  V  cos  <(> 
c  dt  c 


(16) 


This  force  must  be  integrated  over  the  surface  area  of  the 
sphere  which  is  in  contact  with  the  target  to  obtain  the 
total  drag  force, 


Fdrag 


2tt p t  R^sin  <(>  cos^$  d$ 


or 


Fdrag 


sin4 


,1  +  cos 


(17) 


■! 


A 

] 


1 


9 


2  2 

The  quantity  tt  R  sin  <£  is  the  projected  area  in  the 

S 

direction  of  motion  and  is  denoted  A,  .  Hence, 


Fdrag  "  (X  +  cos2<i)s)  A0 


ptVc 


(18) 


Upon  comparing  Eq.  (18)  with  the  usual  definition  for  the  drag 
coefficient 


(19) 


2 

it  is  clear  that  ■  1  +  cos  $  .  For  the  case  of  a  fully 
embedded  sphere,  4>gmv/2  and  therefore  -  1  . 

The  Newtonian  flow  model  is  strictly  valid  only  at  the 
outer  edge  of  the  disturbed  region  and  not  at  the  projectile 
surface.  Hence,  the  value  of  the  drag  coefficient  deduced 
above  must  be  modified  to  account  for  centrifugal  effects  in 
the  distrubed  layer.  It  can  be  shown  that  when  centrifugal 
effects  are  included,  the  pressure  at  the  surface  of  the 
projectile  is  given  by 


P 


pt  Vc  (cos2^  -  |)  . 


(20) 


where  f  is  the  thickness  of  the  layer.  The  layer  thickness 
can  be  estimated  using  a  conservation  of  mass  statement 
together  with  an  assumption  for  the  velocity  profile.  The 
simplest  assumption  is  that  the  velocity  is  uniform  in  the 
layer.  In  this  caBe,  it  can  be  shown  that  S/R  ■  1/2  and, 
further,  that  CD  •  1/2  for  a  fully  submerged  sphere.  This 
is  the  same  value  as  that  obtained  from  the  inertial  term 
in  Eq.  (11)  when  a  -  1  and  4>8  -  tt/ 2  .  The  value  of 
0^  ■  1/2  provides  a  very  good  fit  to  the  data  for  a  wide 

range  of  velocities  and  materials  as  will  be  demonstrated 
shortly.  The  Newtonian  approximation  for  the  drag  term 
probably  works  as  well  as  it  does  because  the  strength  of  the 
target  material  confines  the  region  of  plastic  flow  to  a 
narrow  layer  close  to  the  surface  of  the  projectile  -  a  basic 
assumption  of  the  Newtonian  model. 
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3.  Elastic  energy  per  unit  mass  -  E*e  -  In  Ref.  1, 

the  term  E*  does  not  appear  in  the  momentum  equation. 
Instead,  theelastic  energy  was  included  in  the  model  as  a 
velocity  cutoff  during  the  integration.  This  velocity  cutoff 
was  denoted  by  V*  ,  and  it  was  assumed  that  at  this  velocity 
all  the  remaining  kinetic  energy  of  the  projectile  could  be 
elastically  absorbed  by  the  target  with  no  further  permanent 
penetration.  A  disadvantage  of  using  this  approach  is  that 
Vw  depends  not  only  on  the  physical  properties  of  the  target 
but  also  on  the  density  of  the  projectile.  The  quantity  E*e  , 

however,  depends  only  on  the  physical  properties  of  the  target 
material  ana  is  therefore  a  more  convenient  parameter  to  use. 

During  the  course  of  this  program,  experimental  and 
theoretical  evidence  has  been  developed  to  show  that  E*e 

depends  on  the  dimensionless  penetration  depth  z/2R  . 

However,  in  analyzing  experimental  data,  only  integral 
quantities  such  as  projectile  mass,  radius,  velocity,  kinetic 
energy,  and  final  penetration  depth  are  available.  Empirical 
information  about  the  time  dependence  of  velocity  and  pene- 
tration  are  generally  not  available.  Hence,  the  value  of 
E*e  ,  which  Is  a  local  property,  can  only  be  deduced  from 

impact  data  in  an  integral  sense.  The  method  by  which  this 
can  be  done  is  as  follows. 


Consider  the  equation  for  the  rate  at  which  work  is  done 
on  the  target  (Eq.  (5)).  This  equation  can  be  integrated  to 
give  the  total  work  done  on  the  target 


Adz 


(21) 


where  p  is  the  final  penetration.  The  integral  can  be 
written  as  follows 


W 


t 


where 


E*e<P) 


dz  +  pfc 


1 

v(p) 


E*e(z)Adz  . 


v(p) 


(22) 


(23) 
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The  quantity  E*e  (p)  represents  the  average  value  of  E*g 

over  the  trajectory  and  only  depends  on  the  final  value  of 
the  penetration  p  ,  It  will  be  shown  shortly  that  E* 

is  constant.  Hence  the  quantity  E*p  4-  E,Ve  in  Eq .  (22)  is 

solely  a  function  of  p  .  Therefore .  for  each  impact  test 
a  unique  value  of  E *  -  E^  +  E*e  can  be  found  wnich  when 

integrated  via  Eqs.  (13)  and  (14)  gives  the  correct  pene¬ 
tration  depth  p  for  the  given  impact  velocity  V0q  . 

It  will  be  shown  that  E*(p)  deduced  in  this  manner  can 
be  interpreted  as  the  superposition  of  two  components;  a 
constant  component,  E*  ,  for  relatively  large  penetrations 

and  a  depth-dependent  component,  E*e>  for  shallow  penetra¬ 
tions.  It  was  shown  in  Ref.  1.  that  E*e'  has  the  following 
form 

-  k  <p/2R)N  ,  (24) 

where  k  and  N  are  constants.  The  exponent  N  is 
approximately  -.75  for  most  metals.  For  other  materials, 

N  is  between  -.75  and  -1.5  . 

To  summarize,  Impact  data  (crater  depth  measurements)  are 
used  to  deduce  E*(p)  .  The  constant  portion  of  E*  obtained 
for  large  penetrations  is  E*  .  For  shallow  penetrations, 

Eq.  (24)  is  used  to  evaluate  Che  constants  k  and  N  which 
best  match  the  data.  In  this  way  E>Ve  is  obtained. 

Note  that  this  procedure  requires  knowledge  of  the  final 
penetration.  In  the  general  case  when  p  is  not  known 
a  priori,  the  local  value  of  E*e  and  not  the  average  value 

must  be  used.  The  local  value  can  be  obtained  by  differentia¬ 
ting  Eq.  (23) 

**.<»>  -  5^<p>  +  $$  3?  [Wp>]  •  C25) 

This  value  can  then  be  used  in  the  integration  of  the  momentum 
equation.  (note  the  bar) ,  because  it  is  constant  over 

the  trajectory  and  depends  only  on  the  final  penetration  p  , 
is  considerably  more  convenient  to  use  than  the  local  value 
of  E*e  .  It  will  be  used  exclusively  in  this  report.  The 

bar  will  be  dropped  from  the  notation,  and  the  local  value 
of  E*e  will  no  longer  be  used.  This  approach  permits  the 
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integration  of  Eq .  (10)  as  was  noted  earlier. 

B,  Experimental  Program 

This  section  describes  the  experimental  program  which 
was  used  to  evaluate  the  impact  properties  E*p  and  E^e 

of  the  various  target  materials.  It  contains  a  brief 
description  of  the  A.R.A.P.  Impact  Facility,  a  summary  of 
the  test  matrix,  and  an  illustration  of  the  method  by  which 
E,Vp  and  E*g  are  evaluated  from  impact  dota. 

1.  A.R.A.P.  impact  facility  -  Figure  2  depicts 
schematically  the  A.R.A.P.  Impact  Facility.  This  facility 
consists  of  a  mounted  rifle  and  enclosed  test  tube  and  test 
chamber.  Two  guns*  were  used  for  this  program;  a  Winchester 
.270  caliber  smooth  bore  rifle  and  a  Power  Line  880  Air  Gun. 

The  Winchester  is  permanently  mounted  to  a  fixed  support 
and  is  bore-sighted  on  the  target.  Cartridges  are  hand 
loaded  using  Hercules  2400  gunpowder.  The  projectiles  are 
0.250-inch  diameter  balls  made  of  either  tungsten  carbide 
or  chrome  steel  (AISI  52100) .  The  ball  is  fired  using  a 
Lexan  sabot  which  is  manufactured  in  four  sections.  The 
sabot  separates  aerodynamically  after  leaving  the  muzzle  and 
is  stripped  from  the  flight  path  ahead  of  the  test  chamber. 

The  velocity  range  for  steel  and  tungsten  carbide  projectiles 
is  700  to  5,500  feet  per  second,  Lighter  projectiles  have 
been  fired  at  velocities  up  to  6,600  feet  per  second, 

The  Power  Line  gun  is  utilized  to  extend  the  low  velocity 
range  of  the  facility,  The  gun  can  fire  0.156  and  0.172-inch 
diameter  tungsten  carbide  ana  chrome  steel  balls  in  a  velocity 
range  between  160  and  740  feet  per  second, 

Projectile  velocity  is  measured  using  a  Schmidt-Weaton 
Chronograph.  Two  light  screens,  separated  by  a  fixed  distance, 
sense  the  passage  of  the  projectile  UBing  photo-resistor 
elements.  The  flight  patn  is  illuminated  by  light  which 
passes  through  slits  in  the  shelf  located  in  the  chamber. 

The  shadow  produced  by  the  projectile  on  the  light  screens 
triggers  and  then  stops  a  counter.  Digital  readout  of  the 
velocity  is  provided  on  a  display  board. 

Targets  are  mounted  in  a  cylindrical  tube  at  the  aft  end 
of  the  test  chamber.  MoBt  of  trie  targets  were  circular  disks 


The  Impact  Facility  also  includes  both  a  Weatherby  .460  rifle 
and  a  30-06  rifle,  although  these  rifles  were  not  used  for 
this  test  program. 
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with  a  nominal  diameter  of  6  inches  and  a  thickness  of  1  inch 
Other  shapes  and  sizes,  however,  can  be  easily  accommodated. 
The  target  is  restrained  on  the  front  by  an  annular  steel 
ring  fastened  to  the  cylindrical  steel  tube  and  on  the  rear 
by  a  1-inch  think  wood  disk  and  several  steel  disks  held  in 
place  by  a  1-inch  steel  screw. * 

Documentation  for  each  test  is  available  and  a  summary 
of  the  data  is  contained  in  Appendix  A.  Pre-test  measure¬ 
ments  include  thickness,  diameter  and  mass  of  target,  mass 
and  diameter  of  the  projectile  and  photographs  of  the  undam¬ 
aged  target.  Post-test  measurements  Include  target  crater 
depth  -  using  a  dial  indicator,  and  diameter,  target  mass, 
and  projectile  mass  and  diameter.  Photographs  were  also 
taken  of  the  damaged  target, 

2,  Test  program  -  A  total  of  181  impact  tests  were  con¬ 
duct  e3™on~mrHerent  target  materials.  These  materials 
include  seven  metals,  four  alloys,  two  plastics,  two  ceramics 
and  one  composite,  Table  1  provides  a  description  of  these 
materials  and  a  summary  of  the  results.  A  tabulation  of  the 
data  is  contained  in  Appendix  A. 

The  limits  for  the  velocity  range  over  which  each 
material  was  tested  are  based  on  the  following  criteria.  The 
lower  limit  was  set  by  the  chronograph  sensitivity  which  is 
approximately  150  feet  per  second.  The  upper  limit  of  the 
velocity  range  was  set  to  limit  maximum  crater  depth  to 
approximately  half  the  target  thickness  in  order  to  minimize 
backface  effects. 

Note  that  for  most  of  these  materials,  a  free  fall  test 
was  also  conducted  by  dropping  a  0.50- inch  diameter  chrome 
steel  ball  from  the  ceiling  of  the  laboratory  onto  the  target 
From  a  height  of  approximately  14  feet,  the  impact  velocity 
was  computed  to  be  28  feet  per  second. 

3.  Evaluation  of  and  E^e  -  For  each  impact  test, 

a  value  of  E*  was  computed  using  Eqs.  (13)  and  (14)  which 

fave  the  correct  value  of  penetration  depth,  p  ,  for  the 
mpact  velocity.  Figure  3  shows  the  results  of  these  compu¬ 
tations  for  the  soft  aluminum  target.  The  computed  E,v  is 
shown  as  a  function  of  the  normalized  penetration  p/d  , 


9T" 

The  breakwires  shown  in  Figure  2  are  used  to  measure  projec¬ 
tile  residual  velocity  for  finite  thickness  targets.  The 
wires  were  not  used  for  this  test  program. 
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where  d  is  the  initial  diameter  of  the  ball.  For  this 
material;  the  values  of  E,v  fall  within  a  narrow  band,  and 
the  average  value  of  the  band  is  approximately  84  Btu/lbm. 

All  of  the  E,v  values  fall  within  15%  of  the  average,  except 
for  one  shallow  penetration  data  point  which  was  relatively 
inaccurate. 

As  noted  earlier,  for  each  material  an  equation  of  the 

form 

E*  -  E*p  +  E*g  -  E*p  +  k  (p/dQ) N  ,  (26) 


was  fit  to  the  E*  computations.  The  functional  form  of 
E*  is  based  on  the  static  experiments  described  in  Ref.  1. 
Therefore,  for  soft  aluminum,  E^p  ■  84  Btu/lbm  and  E*e  ■  0  , 

Appendix  B  contains  the  E*  evaluations  for  the  remaining 
materials , 

C.  Data  Theory  Correlations 

After  E,Vp  and  E#e  have  been  evaluated  for  a  given 

target  material,  the  Integral  Theory  of  Impact  can  be  used 
to  compute  the  depth  of  penetration  for  any  impact  velocity. 
The  procedure  is  simply  to  differentiate  E^e(p)  ,  evaluate 
E*e(z)  according  to  Eq.  (25)  and  integrate  Eq.  (13)  and 
Eq.  (14)  simultaneously.  In  this  section,  the  impact  data 
for  each  material  are  presented  and  cormared  to  tne  results 
of  such  computations.  Data  are  shown  for  both  tungsten  car¬ 
bide  and  chrome  alloy  steel  (AISI-52100)  balls, 

In  each  of  the  tests,  except  those  noted  "ball  broken", 
there  was  no  measurable  plastic  deformation  of  the  ball 
after  impact.  For  most  of  the  tests,  the  ball  rebounded 
from  the  target.  Hence,  the  crater  depth  could  be  measured 
directly  using  a  dial  indicator,  In  some  tests,  however, 
the  ball  remained  embedded  in  the  target.  For  these  cases, 
the  dial  indicator  was  used  to  locate  the  position  of  the 
top  of  the  ball  and  the  ball  diameter  was  added  to  obtain 
crater  depth. 

1,  Metals  -  The  results  for  six  metal  targets  are  shown 
in  Figures  4-9.  Each  of  the  targets  was  a  circular  disk 
with  a  nominal  diameter  of  six  inches.  The  nominal  thickness 
of  each  target  was  one  inch,  except  for  a  few  of  the  lead 
targets  which  were  two- inches  thick.  With  the  exception  of 
the  cast  iron  target,  these  metals  are  relatively  soft 
(E*p  <  100  Btu/lbm)  and  inelastic  (E*g  ■  0)  . 
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ALUMINUM  TARGET 
(1100-  F  PLATE) 
E#p  =  84  Btu/lbm 
E  #  e  8  0 


DATA 

O  .250"  WC  BALL 
O  .156  "  WC  BALL 
C  .156"  STEEL  BALL 


A 

4 


.500" STEEL  BALL 
BALL  EMBEDDED 


i_J _ i _ I _ l_J 

100 

Velocity,  ft/sec 


Figure  4 


CADMIUM  TARGET 
99.9  %  PURE,  OPEN  CAST 
E#  p  -  29  Btu/lbm 
£#•  B  0 


10  100  1000 

Velocity,  ft/sec 

Figure  5 
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COPPER  TARGET 
H.R.  ETP  PLATE 
E*p  *38  Btu/lbm 
E*e  =0 

O  .tso"  wc  ball 

O  .156"  WC  BALL 

A  .500"  STEEL  BALL 

□  .156"  STEEL  BALL 

&  BALL  EMBEDDED  IN 
TABOET 

•  P/f  >  0.5  / 


^  WC  bai I  theory 
Steel  ball  theory 


Velocity,  ft/sec 


Figure  6 


IRON  TARGET 

GRAY  CAST  IRON -CLASS  40 
E#p»  12  I  Bfu/lbm 


=  7  (p/d0) 
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O  .»80“  WC  BALL 

O  .186“  WC  BALL 

A  .800“  8TEEL  BALL 

O  .186“  STEEL  BALL 

&  BALL  EMBEDDED  IN 
TAHOE T 

•  BALL  BROKEN 


WC  ball  theory 


Steel  ball  theory 


Velocity,  ft/sec 


Figure  7 


ZINC  TARGET 
99.9  %  PURE,  OPEN  CAST 
E*p  =  62  Btu/lbm 

E  #t*0 


'il  -i - 1 - 1  i- 1  - 

IOO  1 000 

Velocity,  ft/sec 


Figure  9 


In  addition  to  these  metals,  four  tests  were  conducted 
on  silicon  targets.  The  silicon  was  in  the  form  of  relatively 
large  chunks  -  major  dimension  approximately  four  to  five 
incnes  -  obtained  as  a  byproduct  of  steel  processing,  The 
samples  were  porous  and  contained  a  small  amount  of  iron 
impurities.  In  order  to  test  this  material,  each  of  the 
four  largest  pieces  was  cast  in  concrete,  and  the  front 
surface  of  each  piece  was  ground  flat.  Suitable  material 
was  available  for  only  four  tests.  The  average  value  of  E* 
for  these  tests  was  134  Btu/lbm.  The  data  ware  insufficient 
to  determine  E*  . 

2,  Metal  alloys  -  Pour  metal  alloys  were  tested,  and  the 
results  are  shown  in  Figures  10-13.  The  aluminum  alloy  was 
designated  5083  Aluminum,  but  its  temper  was  unknown.  However, 
it  had  a  static  hardness  of  Rockwell  B-75.  The  targets  had 
a  rectangular  cross-aaction,  3.9"  by  4.5",  and  were  one-inch 
thick.  This  alloy  has  a  relatively  high  value  of  E* 

•  245  Btu/lbm  and  also  has  a  moderate  value  of  E*  p  for 
small  penetrations.  e 

Two  steels  were  tested  -  a  mild  steel  (1020  -  hot  rolled) 
and  an  armor  steal  (rolled  homogeneous  steel  armor) .  The  mild 
steel  targets  were  six- inch  diameter  disks i  the  armor  steel 
targets  were  3.9"  by  4,5"  rectangles.  All  steel  targets  were 
one- inch  thick.  The  value  of  E>Vp  varied  from  141  Btu/lbm 

for  the  mild  steel  to  203  Btu/lbm  for  the  armor  steel.  Both 
steels  have  moderate  E*e  components  with  the  armor  steel 

having  approximately  twice  the  value  of  the  mild  steel.  Note 
that  the  tungsten  carbide  balls  were  broken  for  Impact  veloc¬ 
ities  in  excess  of  2,000  feet  per  second,  However,  the 
correlations  with  nondeforming  ball  theory  are  still  fairly 

5ood.  This  is  because  the  tungsten  carbide  exhibits  very 
ittle  deformation  prior  to  brittle  fracture  and  not  more 
than  three  or  four  pieces  are  formed  for  velocities  of  2,000 
feet  per  second.  However,  when  the  ball  is  shattered  into 
many  small  pieces,  the  rigid  ball  theory  overpredicts  the 
penetration  as  is  evident  in  Figure  11  for  velocities  in 
excess  of  4,000  feot  par  second. 

The  titanium  alloy  (Figure  13)  was  T1-6AA-4V.  The 
targets  had  a  square  cross  section  4.375"  on  a  side  and  a 
nominal  thickness  of  0.68  inches.  This  alloy  has  the  highest 
value  of  E*p  of  any  of  the  metals  and  alloys  tested.  It 

also  has  a  moderate  value  of  E*e  . 
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3.  Plastics  -  Two  plastics  -  acrylic  and  polycarbonate  - 
were  tested,  ancT the  results  are  shown  in  Figures  14  and  15. 
Targets  for  each  material  were  six-inch  diameter  disks  and 
one- inch  thick.  The  acrylic  was  Acrylite  and  is  produced  by 
the  American  Cyanamid  Company.  The  polycarbonate  was  Lexan , 
produced  by  the  General  Electric  Company.  These  materials 
have  similar  impact  properties.  The  values  of  E*p  are  not 

high  (~  100  Btu/lbm)  and  differ  by  only  about  10%.  They 
each  exhibit  large  values  for  E*  .  For  penetrations  less 
than  one  ball  diameter,  the  E*g  “component  is  larger  than 

.  Despite  the  similarities  in  the  impact  properties, 

the  appearance  of  the  damaged  targets  is  quite  different. 

The  acrylic  is  brittle  and  exhibits  large  cracks  over  most 
of  its  surface  as  well  as  backface  spall,  The  polycarbonate 
is  much  more  ductile,  and  the  damage  is  localized  to  a  small 
region  in  the  vicinity  of  the  ball. 

4.  Ceramics  -  Two  ceramics  -  glass  and  salt  -  were  tested, 
and  the  results  are  shown  in  Figures  16  and  17.  The  glass  targets 
were  Corning  Pyrex  7740  cast  mirror  blanks,  six  inches  in  diameter 
and  one-inch  thick.  This  material  is  extremely  brittle  and,  as 

a  result,  there  is  considerable  scatter  in  the  crater  depth 
measurements . 

The  salt  targets  which  were  cut  from  large  blocks  of 
polycrystalline  sodium  chloride  were  six-inch  diameter,  one- 
inch  thick  disks.  This  material  is  also  quite  brittle  and 
exhibits  considerable  scatter  in  the  data.  However,  a  value  of 
E*  -  86  Btu/lbm  appears  to  provide  good  correlation  with 
the  data. 

5.  Composite  -  Figure  18  shows  the  results  for  the 
Kevlar  target.  TKis  composite  is  a  woven  fabric  produced  by 
DuPont  and  consists  of  aramid  (Kevlar)  fibers  treated  with 
an  epoxy  resin  and  molded  into  a  rigid  sheet,  The  sheets 
which  are  approximately  3/8"  thick  were  cut  into  four-inch 
squares.  Targets  were  made  by  stacking  three  pieces.  No 
bonding  between  layers  was  employed.  Kevlar  has  a  value  of 
E*p  comparable  to  steel  and  a  large  component  of  E*fi  . 

D.  Summary  of  Qualification  Tests 


A  summary  of  the  value  of  E*  and  the  constants  needed  to 

evaluate  Ew@  is  contained  in  Table  1.  In  general,  the  Integral 

Theory  of  Impact  when  used  with  these  two  properties  provides 
very  good  correlation  with  the  data  for  the  wide  range  of  materials 
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Figure  14 
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which  have  just  been  discussed.  The  theory  correctly  predicts 
the  dependence  of  penetration  depth  on  velocity,  projectile 
diameter,  and  density.  Hence,  the  procedure  described  in 
this  chapter  can  be  used  to  obtain  the  impact  properties  of 
most  target  materials. 

It  appears,  however,  that  the  procedure  does  require  some 
modification  for  characterizing  extrememly  hard,  brittle 
materials.  When  the  target  is  hard,  the  projectile  breaks  at 
relatively  low  velocity  ,  and  insufficient  rigid  ball  data  can 
be  obtained  to  deduce  E*  .  If  only  broken  ball  data  are 
available,  then  one  must  be  careful  to  separate  the  target 
E*  from  the  ball  deformation  model.  Even  with  nondeforming 
balls,  brittle  targets  Inherently  have  large  scatter  in  crater 
depth  measurements,  A  quantitative  measure  of  this  scatter 
or  a  procedure  to  reduce  the  scatter  is  necessary  if  much 
data  for  brittle  materials  are  needed.  Despite  these  limita¬ 
tions  ,  the  procedure  described  in  this  chapter  does  provide 
an  adequate  dynamic  screening  test  for  target  materials  in 
impact  applications. 
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III.  IMPACT  PROPERTIES  OF  MATERIALS 

This  chapter  contains  a  summary  of  the  formulas  which 
relate  the  characteristic  properties  E*p  and  E#e  to  more 

fundamental  material  properties  such  aB  Brinell  hardness, 
melting  temperature,  heat  capacity,  and  Young's  modulus.  The 
derivation  of  these  formulas  was  presented  in  Ref,  1,  and 
will  not  be  repeated  here.  The  values  of  the  constants  which 
appear  in  these  formulas  differ  slightly  from  those  which 
were  contained  in  Ref.  1.  The  revisions  are  the  result  of 
a  more  detailed  survey  of  materials  data  and  a  better  fit 
to  more  impact  test  data.  Also  Included  in  this  chapter  is 
a  comprehensive  table  which  summarizes  the  pertinent  pro¬ 
perties  for  many  target  materials. 

A.  Formulas  for  E*p  and  E*e 

Before  writing  the  equations,  it  is  particularly  impor¬ 
tant  that  the  caveats  be  mentioned.  When  simple  formulas  are 
presented  which  purport  to  give  solutions  to  complex  problems 
the  temptation  is  to  accept  the  results  without  hesitation 
and  to  Ignore  the  limitations  of  the  model  which  produced 
the  formulas.  The  formulas  which  are  written  below  were 
derived  to  serve  as  a  tool  for  the  preliminary  screening  of 
target  materials  for  particular  impact  applications,  not  as 
a  final  solution  for  and  E*e  and  not  as  a  substitute 

for  impact  tests. 

If  one  follows  the  derivations  in  Ref.  1.  closely,  it 
becomes  clear  that  some  grosB  assumptions  have  been  made 
regarding  material  behavior  and  failure  mechanisms,  and  that 
oraer-of -magnitude  arguments  have  been  used  to  justify  the 
use  of  a  single  equation  which  is  applicable  to  hard  steel  as 
well  as  soft  zinc  and  to  brittle  ceramic  as  well  as  ductile 
plastic.  The  hope  was  that  such  a  model  would  give  a  ball¬ 
park  estimate  for  E^p  and  E*e  .  In  this  regard,  the 

formulas  succeed  as  a  tool  to  distinguish  good  target  mater¬ 
ials  from  poor  materials  and  to  provide  a  rough  ranking  of 
materials.  The  fact  that  the  equations  predict  the  values 
of  E*p  and  E*fi  within  157.  for  many  materials  is  a  dividend 

and  a  confirmation  that  the  assumptions  of  the  model  are  good 
for  many  materials.  More  importantly,  however,  the  equations 
isolate  those  material  properties  which  have  the  most 
influence  on  impact  performance. 
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1.  Ey,e  -  The  formula  for  E*fi  Is  given  by 


£ 


*e 


.69  x  10"3 


(27) 


o 

where.  B  Is  Brinell  hardness  in  N/m  ,  p  is  maBs  density 

in  kg/m  ,  E  is  Young's  modulus  in  N/m  ,  y  is  a  dimen¬ 
sionless  strain  rate  parameter,  and  E*e  is  given  in 

Btu/lbm.  The  parameter  y  accounts  for  the  increase  in 
Brinell  hardness  at  high  strain  rates.  It  can  be  evaluated 
from  hardness  measurements  taken  at  room  temperature  and  at 
crvogenic  temperatures  using  the  strain  rate  -  temperature 
relationships  discussed  in  Ref.  1  and  in  Chapter  13  of 
Ref.  5.  Typically,  y  has  a  value  of  1.5  for  metals  and 
5.0  for  plastics . 


2. 


-  The  equation  for 


C  T 

_g _ m 

2326. 


Jin 


E*p  is  given  by 

+  j] 


(28) 


where  Cp  is  the  specific  heat  in  joule/kgK,  Tm  is  the 

melting  temperature  in  K  ,  af  is  the  flow  stress  corrected 

for  strain  rate  effects  in  N/m2  ,  E*p  is  in  Btu/lbm,  and  C2 

and  8  are  constants  which  are  obtained  by  correlating  this 
equation  with  data  from  impact  experiments.  T^  is  the  initial 

temperature  of  the  material,  and  e  1b  the  strain  rate.  The 
present  best  estimate  for  the  constants  is 


S  -  .072 

C2  -  5.76  .  (29) 

Note,  the  factor  2326  in  the  denominator  of  Eq.  (28)  Is 
required  to  convert  units  from  joule/kg  to  Btu/lbm  -  the 
typical  units  in  which  E*p  is  expressed.  The  quasi- 

static  flow  stress  can  be  obtained  from  Brinell  hardness 
measurements  using  the  Frandtl  solution  (Ref.  5)  for  the 
deformation  flow  field  in  a  hardness  test  given  by 


5.  HcClintock,  F.  A.  and  Argon,  A.  S.,  Mechanical  Behavior  d 
Matciiali ,  Chapter  13,  Addison-Wesley ,  Reading,  Mass., 

1966 
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(30) 


This  flow  stress  is  corrected  for  strain  rate  effects  using 
the  "velocity-temperature"  or  "temperature-strain  rate"  inter¬ 
relationship  which  was  mentioned  in  Ref.  1  and  is  described  in 
Ref.  5.  In  essence,  this  correction  eauates  the  hardness  test 

Performed  at  high  strain  rate  to  a  hardness  test  performed  at 
over  temperature.  For  many  materials  the  Increase  in  hardness 
for  impact  strain  rates  is  equivalent  to  the  increase  in  hard¬ 
ness  for  a  100  kelvin  decrease  in  temperature.  In  general, 

Eq.  (28)  is  within  approximately  15%  of  the  experimentally 
deduced  values  of  for  most  metals. 

For  the  brittle  materials,  the  theory  is  generally  within 
a  factor  of  two  of  the  experimental  data.  The  discrepancy 
between  theory  and  data  is  not  unexpected.  Recall  that  EWp 

is  the  plastic  energy  dissipation  which  occurs  in  a  thin, 
high  shear  region  adjacent  to  the  projectile.  In  impact 
tests  of  brittle  targets,  there  is  a  large  shard  fracture 
mode  which  limits  the  extent  of  the  shear  region.  In  essence, 
the  tensile  release  waves  at  the  side  of  the  projectile  break 
off  large  pieces  of  target  material  before  significant  shearing 
can  occur,  thereby  inhibiting  the  energy-dissipation  potential 
of  the  material.  Hence  the  performance  of  brittle  materials 
in  impact  tests  is  considerably  less  than  the  theoretical 
prediction. 

Equation  (28)  has  one  other  interesting  feature  which 
should  be  noted.  When  the  value  of  o^/gp  Tm  is  very 

small  compared  to  1,  Eq.  (28)  can  be  written 


(31) 


For  this  case,  there  is  a  linear  relationship  between 
E*p  and  the  Brinell  hardness  of  the  material.  However,  when 

the  value  of  a^/gp  Cp  Tm  is  large,  then  the  logarithm  cannot 

be  simplified  and  increasing  the  hardness  of  the  material  may 
result  in  but  small  improvements  in  E*p  because  of  the 

nature  of  the  logarithm  function,  Hence,  to  a  first  approxi¬ 
mation,  Eq.  (28)  can  be  used  to  determine  those  materials 
whose  impact  performance  can  be  improved  by  hardening  techniques. 
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B.  Table  of  Materials  Properties 

Table  2  presents  a  listing  of  the  pertinent  materials 
properties  of  many  materials  for  which  hardness,  specific 
heat,  and  melting  temperature  data  are  available.  These 
data  were  obtained  from  a  variety  of  sources  in  the  open 
literature.  The  properties  were  used  in  Eq.  (28)  to  obtain 
the  value  of  E*,.  . 

The  final  column  in  Table  2  lists  the  adiabatic  hardness  - 
the  product  pE*p  -  of  each  material.  This  quantity  is 

particularly  Important  for  determining  the  penetration  of  a 
deforming  projectile.  Recall  from  Eq,  (9)  that  the  pressure 
on  the  front  race  of  a  projectile  is  primarily  pE*p  for 

low  velocities.  This  pressure  should  be  as  high  as  possible 
if  the  objective  is  the  deformation  of  the  projectile. 

Note  that  the  properties  in  Table  2  contain  both  metric 
and  English  standard  units.  Both  E*p  and  the  product  pE#p 

have  traditionally  been  tabulated  in  the  units  shown  in  Table 
2.  To  convert  these  quantities  to  their  metric  equivalent  the 
following  conversion  factors  may  be  employed. 
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IV.  LONG  ROD  PENETRATOR  MODEL 


This  chapter  summarizes  the  integral  theory  for  long 
rod  penetrators.  Preliminary  development  of  this  model  was 
described  in  Ref.  2,  and  listings  of  the  computer  codes  ROD 
and  PEN  were  provided.  Unfortunately,  these  codes  contained 
some  theoretical  and  numerical  errors  which  produced  anomalous 
results  at  times.  Improvements  have  been  made  to  the  ROD 
code  which  eliminate  much  of  the  anomalous  behavior.  However, 
it  is  important  to  note  that  there  are  still  certain  limita¬ 
tions  in  the  use  of  this  code.  In  its  present  form,  the  model 
is  limited  to:  (1)  high  velocity  impact,  (2)  homogeneous, 
isotropic  rods,  and  (3;  rod  materials  with  some  ductility. 

The  reasons  for  these  limitations  are  described  below. 

In  the  model,  the  head  of  the  rod  is  assumed  to  flow 
hydrodynamlcally.  There  are  no  constitutive  equations  built 
into  the  deformation  model,  and  there  is  no  mechanism  for 
energy  dissipation  due  to  internal  shear  stresses.  As  a 
result,  the  model  is  applicable  for  Impact  velocities  which 
are  high  enough  so  that  the  pressures  generated  by  the  Impact 
are  much  In  excess  of  the  material  yield  strength. 

An  additional  assumption  of  the  model  is  that  the  shape 
of  the  flowfield  in  the  head  does  not  change  very  much  from 
material  to  material.  Hence,  each  material  undergoes  the 
same  amount  of  plastic  deformation  prior  to  failure.  The 
model  parameter,  e0  ,  which  represents  the  plastic  strain  in 

the  head,  accounts  for  this  deformation.  A  single  value  for 
eQ  appears  to  be  adequate  for  most  ductile  rod  materials. 

However,  for  brittle  materials  which  have  much  lower  values 
for  strain  to  failure,  a  different  model  for  r  may  be 
needed. 

A  further  assumption  of  the  model  is  that  the  rod  material 
is  homogeneous  and  isotropic.  Anisotropic  material  properties 
such  as  might  exist  for  composite  rods,  for  example,  cannot 
be  handled  In  the  present  model.  In  the  future,  as  more 
Impact  data  become  available,  the  model  will  be  extended  to 
Include  constitutive  equations  and  anisotropic  material 
properties. 

In  what  follows,  the  equations  which  govern  the  penetra¬ 
tion  of  a  long  rod  projectile  will  be  presented.  These  equa¬ 
tions  form  a  coupled  set  of  ordinary  differential  equations 
and  are  solved  numerically  using  the  ROD  code.  Appendix  C 
contains  the  computer  listings  and  a  user's  guide  for  this 
code.  Typical  results  of  numerical  computations  are  pre¬ 
sented  at  the  end  of  this  chapter, 
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A.  Physics  of  the  Model 


i 


i 

I* 
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It  is  known  from  x-ray  photographs  that  the  stages  of 
long  rod  penetration  may  be  roughly  characterized  as  in 
Figure  19.  As  the  rod  impacts  tne  target,  the  pressures 
generated  at  the  interface  begin  to  deform  the  front  end  of 
the  rod,  as  in  Figure  19(b).  Simultaneously,  the  same 
pressure  erodes  the  target  and  produces  a  crater  as  in 
Figure  19(c).  As  penetration  continues,  material  at  the 
leading  face  of  the  projectile  is  eroded  by  the  target, 
forced  out  laterally  from  the  contact  region  by  the  high 
pressure  there  and  ejected  back  out  of  the  crater.  As 
material  erodes  from  the  rod  face,  new  material  is  supplied 
to  this  region  by  the  shaft  of  the  rod  which  is  traveling 
at  a  higher  velocity  than  the  rod-target  interface.  At 
some  point,  the  shaft  material  is  completely  used,  Figure  19(d), 
and  the  remaining  head  is  then  brought  to  a  stop  by  the 
target. 

1,  Kinematics  -  The  model  for  the  rod  flowfield  consists 
of  two  regions;  the  head  of  the  rod  which  corresponds  to  the 
region  undergoing  hydrodynamic  strain  and  the  shaft  or  rear 
portion  of  tne  rod  which  contains  the  undeformed  material. 

During  penetration  the  head  which  is  in  contact  with  the 
target  decelerates  and  spreads  laterally,  It  is  assumed 
that  the  rod  is  cylindrical,  that  the  shape  of  the  head  remains 
cylindrical  during  the  penetration,  and  that  the  velocity 
flowfield  in  the  head  is  linear.  Figure  20  schematically 
illustrates  the  rod  model  and  the  pertinent  nomenclature. 


The  motion  of  the  material  in  the  head  is  described  by 
a  center  of  mass  velocity  Vcm  and  by  the  velocity  of  its 

front  face  and  side  face  Vt  .  At  any  instant,  the  thick¬ 

ness  of  the  head  is  21  ,  and  its  radius  is  b  .  Since  the 
velocity  field  is  assumed  to  be  linear,  the  following  equations 
can  be  written  for  the  particle  velocity  at  the  shaft-head 
interface,*  , 


V,  + 


2<Vcm  -  V*> 


(32) 


'1 


;  i 


!  ] 
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and  for  the  rate  of  change  of  head  dimensions 


FLOW  FIELD  OF  ROD 
AT  VARIOUS  STAGES  OF  PENETRATION 


TARGET  TARGET 


<d> 


ACTUAL 

DEFORMATION 

Figure  19 
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OUTFLOW 

assumptions: 

I  .  Conservation  of  energy 

2.  Conservation  of  momentum 

3.  Linear  flow  field  In  head 

4.  Continuity  of  mass  flow  across  interface 

5.  Constant  yield  stress  at  interface 

6.  Mass  of  penetrator  erodes  from  head 
when  radius  exceeds  (l+€0)a 

7.  The  model  depends  upon  two  parameters: 


The  yield  stress  a  *pE+tmXPE* 


and 


The  shearing  radius  given  by(l  +  €0)a 

Figure  20 


(33) 


V*  '  Vcm  • 


b  ■  Vp  for  b  <  (1  +  r.o)a  , 
b  ■  0  for  b  >  (1  +  £Q)a  , 


(34a) 

(34b) 


where  a  is  the  initial  radius  of  the  head,  and  the  dot 
above  a  symbol  denotes  the  time  derivative.  The  quantity  e0 

represents  the  allowable  plastic  strain  in  the  head. 

As  penetration  proceeds,  the  head  widens  as  rod  material 
is  forced  to  flow  in  the  lateral  direction.  At  some  distance 
from  the  axis  of  the  rod,  say  (1  +  cQ)a  laterally,  material 

is  detached  from  the  head.  The  dynamics  of  this  material  as 
it  further  interacts  with  the  target  haB  no  subsequent  effect 
on  the  head  or  shaft.  This  assumption  is  justified  because 
the  rod  material  at  this  point  in  the  flow  has  been 
adiabatically  heated  by  plastic  work  to  such  an  extent  that 
its  shear  strength  is  very  low.  Hence,  it  is  able  to 
influence  the  rod  only  through  compressive  or  hydrodynamic 
forces.  However,  the  axial  component  of  the  compressive 
hydrodynamic  force  will  only  be  significant  within  one  or 
two  rod  radii  from  the  central  axis.  Thus,  (1  +  eQ)a  <  2a 

and  really  characterizes  the  turning  radius  of  the  rod 
material  in  the  target  or  the  shape  of  the  flowfield  in  the 
head.  For  simplicity,  it  is  assumed  that  the  shape  of 
the  flowfield  does  not  change  too  much  with  different  materials 
Hence,  eQ  is  the  same  for  all  rod  penetrators,  independent 

of  material.  Its  value  can  be  deduced  from  impact  data. 

More  on  this  point  later . 

2.  Mass  Conservation  -  Conservation  of  mass  across  the 
shaft-head  Interface  imposes  the  following  condition  for  the 
rate  at  which  material  flows  from  the  shaft  to  the  head, 


ma  -  Pp*a"  <VS  -  Vb)  ,  (35) 

where  pp  is  the  penetrator  density,  and  Vg  is  the  shaft 
velocity. 

When  the  radius  of  the  head  reaches  the  value  (1  +  e0)a  , 
it  is  assumed  that  any  further  increase  in  the  radius  simply 


t-f 


results  in  loss  of  rod  material  across  the  boundary.  The 
rate  of  mass  loss  from  the  head  is  given  by 


mb  ■  4TTppbUVt  -  b)  .  (36) 

* 

Before  the  head  reaches  the  cutoff  radius,  V  •  b  and 
■  0  .  After  the  cutoff  radius  is  reached,  d  ■  0  . 

The  rate  of  change  of  mass  in  the  head  can  be  obtained 
from  Eqs.  (35)  and  (36) 

nih  ■  ma  -  mb  “  upp  a2(VB  -  Vb)  -  4bJ!.(Vt  -  b)l  ,  (37) 


and  the  rate  of  change  of  mass  in  the  shaft  is  simply 

®s  ■  -  -  -  va2<V  V  1  (38) 

3,  Momentum  conservation  «  The  pressure  applied  to 
the  front  face  of  the  rod  by  the  target  is  given  by 

*  ■  Pt(E*  +  °D  7)  '  <39) 

This  pressure  acts  across  the  entire  front  face  of  the  rod 
which  is  in  contact  with  the  target.  The  contact  area  is 
2 

Trb  .  Hence,  the  drag  force  becomes 

F  -  1Tb2Pt  +  CD  yj  .  (40) 


The  equation  for  conservation  of  axial  momentum  can  be  written 


Equation  (39)  is  written  for  a  total  E ,,  -  E 


only 


is  used  at  the  present  time. 


*  r 


However , 
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(41) 


A  <“hVcn  +  ».V  * 


where  the  second  term  on  the  right  side  accounts  for  the  flux 
of  axial  momentum  out  the  side  of  the  head.  This  equation 
may  be  separated  Into  two  equations  for  the  head  and  the 
shaft.  The  shaft  will  only  experience  deceleration  forces 
if  it  has  a  nonzero  yield  strength  o  .  In  general,  the  force 
on  the  shaft ,  F  ,  is  given  by 

5 


f-Vcp. 


-  »bip.lE.  +  C. 


-Ven  ' 


F,  -  -  na2o  -  msVs  .  (42) 


When  o  ■  0  ,  bb  is  the  case  for  a  shaped- charge  jet  which 
is  liquid,  the  shaft  velocity  remains  constant  during  the 
entire  penetration. 

An  illustrative  equation  is  obtained  by  combining  Eqs.  (37), 
(38) ,  (41) ,  and  (42)  to  obtain  an  equation  for  the  deceleration 
of  the  head 


The  first  term  on  the  right  is  the  deceleration  of  the  head 
due  to  the  target  pressure.  The  second  term  is  the  accelera¬ 
tion  of  the  head  due  to  the  push  from  behind  by  the  yield 
strength  of  the  shaft.  The  third  term  is  the  net  momentum 
flux  to  the  head  due  to  mass  flux  across  material  surfaces. 

4.  Energy  conservation  -  The  total  kinetic  energy  in  the 
rod  is  given  Sy 


where 


K  -  Ka  +  +  K 

s  cm  r 


Ks  ■  i  vJ  ■ 


is  the  kinetic  energy  in  the  shaft  and 


(44) 

(45) 
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rmwsi*' 


(46) 


Kcm  "  \  mhVcm  * 

Is  the  kinetic  energy  associated  with  the  mass  center  motion 
of  the  head  and 


Kr  -  S  »h  [l  +  <V1  -  V»>2]  • 


is  the  kinetic  energy  in  the  head  relative  to  the  mass  center 
motion,  l.e.,  the  energy  associated  with  the  deformation 
of  the  head,  Equation  (47)  can  be  obtained  by  integrating 
the  linear  velocity  distribution  over  the  entire  volume  of 
the  head,  i.e. 


Kh  ■  Kcm  +  Kr  ‘  i  »p 


ii 


vr  +  ve  +  vz)  RdRdedz  . 


where 


't  ■  E  vt 


Ve  -  0 


v£  vcm  +  f  -  Vcm>  • 


Equations  (44-47)  can  be  differentiated  to  yield  the  followinp 
kinematic  expressions 


Kcm  +  K„  +  Kr  • 

I  mhVcm  +  ”1hVcmVcm 

iv!*  ».v.  ■ 


•  p  "IP  " 

-T  “Vt+(V1  'Vctn>2  +  TT  3VtVt  +  2(V1-Van)ff1-Von).(52) 

0  2  J  L 


50 


The  rate  at  which  work  is  done  on  the  rod  by  the  target 
is  given  by 

U  =  rrb2ptVA^E,v  +  ^V^  .  (53) 

The  heating  rate,  denoted  by  W  ,  or  the  rate  at  which  rod 
material  is  converted  into  the  hydrodynamic  state  is  given 
by 

W  -  .  (54) 

W  represents  the  rate  at  which  energy  is  dissipated  as  mater¬ 
ial  crosses  the  Interface  between  shaft  and  head  and  is  trans¬ 
formed  from  the  solid  to  the  hydrodynamic  state.  The  quantity 
is  the  "adiabatic  yield  strength"  of  the  rod  material. 

More  about  this  quantity  later. 

The  flux  of  kinetic  energy  which  crosses  the  lateral 
boundaries  of  the  head  is  given  by 

K  m  J”b[V?  +V™  +  J  <V>-Vcm>2]  '  <55> 

»  «  *  » 

Conservation  of  energy  requires  that  K+U+W+K^*0  or 

K  -  -nb2PtV14*  +  ^2v5j-»aE*d 

'  7'ib[V?  +  Vcm  +  5  <V*  -Von>J]  '  <56> 

5.  Incompressibility  -  The  final  equation  which  is 
needed  to  specify  the  problem  is  the  statement  of  incompres¬ 
sibility.  At  any  instant,  the  mass  in  the  head  is  given  by 

m^  ■  27rb2Hpp  .  (57) 
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If  the  rod  material  is  incompressible,  then  p  is  constant 
and  Eq.  (57)  can  be  differentiated  to  yield  P  „ 


ffih 


(58) 


6.  System  of  equations  -  The  preceding  equations  form  a 
system  of  first-order  ordinary  differential  equations  which 
can  be  solved  simultaneously  to  yield  solutions,  as  a  function 
of  time,  for  the  velocity  field,  dimensions  of  the  head  and 
shaft,  and  energy  partitioning.  One  additional  equation  is 
needed  to  obtain  solutions  as  a  function  of  depth  of  penetra¬ 
tion,  i. e .  , 


y  “  Vj 


(59) 


The  complete  system  of  equation?  is  summarized  in  Table  3. 

The  numerical  procedure  employed  to  integrate  this  system  is 
discussed  in  Appendix  C.  The  input  parameters  required  for 
running  the  code  are  the  length  and  diameter  of  the  rod,  the 
density,  adiabatic  yield  strength  and  adiabatic  hardness  of 
the  rod  material,  the  density  and  hydrodynamic  mode  energy  of 
the  target,  and  the  plastic  strain  in  the  head  . 

7.  Model  parameters  -  It  is  assumed  that  the  energy 
dissipation  parameter  of  the  rod  material,  ,  is  directly 

proportional  to  the  hydrodynamic  mode  energy,  E,v  ,  i.e. 


E*d  "  X  E* 


(60) 


The  value  of  E*  ,  of  course,  is  deduced  from  impact  tests 
using  the  procedure  described  in  Chapter  II  or,  in  the  absence 
of  impact  data,  estimated  from  Eq.  (28).  The  value  of  E*j 
is  different  from  E*  because  the  flowfield  in  the  rod  isa 
somewhat  different  tnan  in  a  semi-infinite  target  of  the  same 
material.  The  value  of  x  can  be  obtained  by  correlating 
impact  data  and  numerical  computations.  As  a  result  of  such 
correlations,  x  ■  0.42  has  been  obtained.  An  experimental 
program  was  designed  and  conducted  at  A.R.A.F.  to  measure 

at  laboratory  strain  rates.  This  program  is  described 

in  Chapter  V. 
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TABLE  3. 


ROD  MODEL  EQUATIONS 


Mass 


mh  -  op 


(v,-Vb)-m. 


”b  "  4"bll>p(v0 


m. 


v* 


:(wb) 


Momentum 


d  /  V  2  /  CDvJ\  . 

57  (mhVcm  +  insVs)  "  “  "b  pt^  +  — j-bV. 


cm 


mgVg  ■  •  tts  c 


Energy 


K-  KOT  +  Kr  +  Ks 


K  -  i  m,V^_  +  nvV  V 
cm  2  n  cm  n  cm  cm 


kB-  £  mgVg  +  tn.V.V. 


8  8  8 


Kr-  ^  [|  v2  +<Vi-VCB)2]+^  fv^er.-VM 

V2(vvb) 


°D  2 


K  -  -fib  ptVA  ^E*  +  ~  VA 


'*d 


Kinematics 


-K  [v?+vL+i  (V*-V=J2] 


b  -  Vf 


b  <  (l+e0)a 


V.  -V 
*  cm 


b  -  0 

Si 

vb  ■  vi  +  2<VcnfVl> 
y  -  Vj 

Incompressibility 


b  >  (l+e0)a 


The  adiabatic  yield  stress,  c  ,  is  related  to  the 
adiabatic  hardness,  ppE,v  .  of  the  material  by 


a  =  'pE*d  =  XPpE*  '  (61) 

This  equation,  which  is  similar  to  Eq.  (30),  is  analogous  to 
the  relationship  between  the  Brinell  hardness,  B  ,  and  the 
uniaxial  tensile  strength  of  a  material  in  static  tests. 

The  laBt  parameter,  eQ  ,  is  assumed  to  be  constant  for 

all  rod  materials.  The  limitations  of  this  assumption  for 
brittle  materials  and  composite  rods  have  already  been  noted. 
The  value  of  cQ  can  only  be  deduced  by  correlating  data 

and  computations.  As  a  result  of  such  correlations,  the 
value  c0  ■  0.36  has  been  obtained. 

8.  Backface  Effects  -  The  resistance  of  the  target  to 
penetration  by  the  projectile  is  given  by  Eq.  (40).  When 
the  target  is  thick,  the  value  of  E*  obtained  by  the  pro¬ 
cedures  described  in  Chapters  II  and  III  should  be  used. 
However,  to  apply  the  integral  theory  to  targets  of  finite 
thickness,  it  is  necessary  to  adjust  E*  for  backface 
effects. 

Recall  that  the  E*  concept  was  originally  developed 
for  the  flow  of  target  material  around  a  penetrator  in  a 
semi-infinite  target.  The  shear  work  done  on  the  target 
material  in  the  flow  volume  defines  E*  .  When  the  projectile 
has  penetrated  almost  all  the  way  through  the  target,  to 
within  one  or  two  rod  diameters  of  the  backface,  the  target 
material  can  spall  or  simply  bulge  on  the  backside,  rather 
than  flow  around  the  penetrator  hydrodynamically .  Thus, 
each  small  volume  of  target  material  absorbs  less  energy  than 
it  would  in  the  semi-infinite  case.  Thus,  the  effective  E* 
for  the  target  decreases  near  the  backface.  In  Chapter  IV 
of  Ref.  2,  a  simple  model  for  the  decrease  in  E*  near  the 
backface  was  proposed.  The  reader  is  referred  to  that 
development  for  the  details.  Briefly,  for  those  materials 
which  exhibit  backface  effects,  the  value  of  E*  is  held 
constant  until  the  head  of  the  rod  reaches  a  point  which  is 
two  rod  diameters  from  the  backface.  From  this  point,  the 
E*  of  the  target  is  decreased  linearly  with  distance  such 
that  E*  is  zero  when  the  front  of  the  head  reaches  the 
backface . 
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9.  Oblique  Impact  -  The  rod  model  described  above  can 
also  be  used  for  the  case  of  oblique  impact.  In  this  case, 
the  target  thickness  is  equivalent  to  the  slant  length  of 
the  target.  For  oblique  impact  of  finite-thickness  targets, 
the  backface  effect  occurs  when  the  perpendicular  distance 
from  the  front  of  the  head  to  the  backface  of  the  target 

is  less  than  two  rod  diameters.  This  treatment  of  oblique 
impact  does  not  attempt  to  handle  fracture  of  the  rod 
shaft  or  jetting  of  the  rod  front  end  during  impact. 

10.  Multi-layer  targets  -  The  rod  model  has  been  used 
with  some  success  tor  multi-layered  targets.  For  such 
targets,  the  equations  are  solved  for  each  layer.  The  condi¬ 
tions  which  exist  at  the  rear  of  one  layer  simply  become 

the  initial  conditions  for  the  next  layer.  This  procedure 
continues  until  the  penetrator  is  stopped  or  perforates  the 
last  layer. 


B.  Correlations  of  Theory  and  Data 

The  parameters  cQ  and  x  can  be  deduced  from  correla¬ 
tions  with  data.  Experiments  were  conducted  for  A.R.A.P.  at 
the  Ballistics  Research  Lab  (BRL) ,  Aberdeen  Proving  Grounds. 
Several  long  rods  with  a  length- to-diameter  ratio,  L/D,  of  10 
were  fired  into  targets  of  Rolled  Homogeneous  Steel  Armor.  In 
order  to  avoid  backface  effects,  very  thick  targets  -  target 
thickness  greater  than  twice  the  total  rod  penetration  - 
were  used.  The  rods  were  chosen  to  provide  a  variety  of  materials 
and  strengths.  The  materials  included  mild  (1018)  steel, 
tantalum,  lead,  and  Mallory  3000,  a  tungsten  alloy. 

The  data  were  compared  to  calculations  using  the  ROD  code 

to  select  a  best  fit  for  t  ,  which  characterizes  the  maximum 

o 

strain  in  the  head,  and  x  which  relates  the  adiabatic  hardness 
to  the  adiabatic  yield  strength  of  the  rod.  The  value  of  E* 
for  the  armor  steel  target  is  200  Btu/lbm.  For  the  rod 
materials,  the  values  or  E*  are:  Mallory  3000  -  50  Btu/lbm; 
tantalum  -  80  Btu/lbm;  1018  steel  -  140  Btu/lbm;  lead  -  3.5 
Btu/lbm.  Based  on  these  comparisons,  it  was  deduced  that 
eQ  •  .36  and  x  ■  *42  .  A  comparison  of  the  correlations 

between  theory  and  data  is  shown  in  Figures  21  -  24.  Note  that 
the  high  velocity  lead  rod  deformed  upon  exit  from  the  barrel 
and  had  an  irregular  shape  upon  impact.  The  L/D  for  this  rod 
was  somewhere  between  5  and  10. 

The  ROD  code  can  be  used  for  finite  thickness  targets. 

The  values  of  the  parameters  and  x  remain  0.36  and  0.42 

respectively.  However,  the  decrease  in  the  E*  of  the  target 
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Penetration,  inches 


MILD  STEEL  ROD 
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L/D  =10 
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MALLORY  3000  ROD 
INTO  SEMI-INFINITE  RHA 


Penetration,  inches 


LEAD  ROD 

INTO  SEMI-INFINITE  RHA 


Velocity,  kft/sec 

Figure  24 
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in  the  vicinity  of  the  backface  must  be  taken  into  account,  as 
was  noted  above.  Figures  25-28  show  some  comparisons  of 
theory  and  data  for  residual  velocity  and  residual  mass.  The 
theoretical  prediction  for  residual  velocity  is  the  shaft 
velocity  at  the  instant  the  head  of  the  rod  reaches  the 
backface.  The  residual  mass  is  the  sum  of  the  mass  in  the 
shaft  and  in  the  head  at  this  Instant.  The  data  are  taken 
from  a  recent  BRL  report  by  Lambert  (Ref.  6).  Residual  mass 
data  are  of  two  kinds:  (1)  reduced  from  radiographs;  (2)  mea¬ 
sured  from  recovered  fragments.  The  residual  velocity  data 
were  reduced  from  radiographs.  The  projectiles  were  65  gram 
Bearcat  steel  rods  (hardness  Rockwell  0-5-5)  with  various  L/D 
ratios.  These  rods  ware  fired  into  various  thickness  targets 
of  Rolled  Homogeneous  Steel  Armor. 

Figures  25  and  26  show  the  results  for  L/D  ■  5  rods  fired 
into  two  different  finite  thickness  targets.  Figure  27  shows 
the  residual  velocity  for  an  L/D  ■  10  rod.  No  residual  mass 
data  are  available  for  these  tests.  Figure  28  is  for  an 
L/D  *  20  rod,  In  general,  the  agreement  with  data  is  good, 
and  the  model  correctly  predicts  the  trends  with  velocity. 

For  these  cases,  the  model  underpredicts  the  ballistic  limit 
velocity  by  approximately  10  to  15%. 


6.  Lambert,  John  P. :  The  Terminal  Ballistics  of  Certain 

65  Gram  Long  Rod  Penetrators  Impacting  Steel  Armor  Plate, 
kept.  No.  ARBRL-TR02072 ,  May  1978. 


EXPERIMENTS 


This  chapter  contains  a  description  of  an  experimental 
program  conducted  at  A.R.A.P.  to  measure  the  fundamental 

quantity  E*d  .  Recall  that  this  quantity  represents  the 

dissipated  energy  absorbed  by  the  penetrator  material  as  it 
is  transformed  from  the  solid  to  the  hydrodynamic  state.  In 
essence,  Ev<d  is  to  the  penetrator  material  what  E^  is 

to  the  target  material.  An  experiment  was  designed  by  A.R.A.P. 
to  measure  both  the  total  energy  input  to  a  material  specimen 
prior  to  fracture  and  the  residual  nondissipated  energy  of 
the  specimen  (kinetic  energy  in  the  debris  cloud)  after 
fracture.  The  quantity  E*d  can  be  obtained  from  these 
measurements.  ° 

In  the  next  section,  the  experimental  concept  is 
described  in  more  detail.  In  the  following  section,  the 
experiment  including  apparatus,  instrumentation,  and  test 
procedure  is  summarized.  This  is  followed  by  a  discussion 
of  the  results. 

A.  Test  Concept 

The  long  rod  model  which  was  described  in  Chapter  IV 
contains  the  fundamental  parameter  E*d  which  represents  the 

dissipated  energy  absorbed  by  the  penetrator  material  during 
impact.  In  addition,  the  model  contains  a  parameter  o  , 
the  adiabatic  yield  stress ,  which  is  lelatedto  E*d  (o  ■  p  E*d  , 

see  Eq.  (61)).  Computations  of  long  rod  penetration  are 
sensitive  to  a  and,  therefore,  to  the  value  of  E*d  . 

One  method  to  evaluate  E^d  is  to  note  that  in  the  impact 
process  for  like  materials  the  penetrator  and  target  are 
undergoing  similar  pressures  and  deformations.  Therefore,  the 
energy  dissipation  mechanism  in  the  penetrator  material  must 
be  similar  to  that  in  the  target  material.  The  energy  dissi¬ 
pation  in  the  target  is,  of  course,  given  by  E*  .  Since 
the  flow  fields  in  the  target  and  penetrator  are  somewhat 
different,  one  would  not  expect  that  E*d  ■  E,v  .  However,  it 

is  reasonable  to  postulate  that  E*d  night  be  proportional 

to  E>v  .  The  constant  of  proportionality,  which  is  denoted 
by  x  in  Eq.  (60),  then  becomes  a  model  coefficient  and  ia 
evaluated  by  correlation  of  the  model  with  data.  This  is 
the  procedure  which  was  described  in  the  preceding  chapter. 


Another  method  is  to  attempt  to  measure  E*d  directly. 

An  experiment  was  designed  by  A.R.A.P.  to  measure  each  of  the 
components  of  the  energy  partition  associated  with  penetra- 
tor  deformation  and  fracture.  By  writing  an  energy  balance 
for  the  process,  it  is  possible  to  deduce  the  value  of  E*d  . 

During  the  impact  process,  material  in  the  penetrator 
is  being  rapidly  decelerated  in  the  forward  direction  and 
accelerated  in  the  lateral  direction,  so  an  inertial  term 
which  takes  into  account  the  kinetic  energy,  similar  to  the 

CjjV^ (2  term  in  the  target  material,  will  be  present  in  the 

penetrator.  Also,  the  penetrator  absorbs  elastic  and  plastic 
energy  as  it  deforms.  Hence,  the  energy  balance  prior  to 
fracture  may  be  written 


W  “  K  +  E  +  P  (prior  to  fracture)  ,  (62) 


where  W  denotes  the  total  work  done  on  the  penetrator  by 
the  target,  K  is  the  kinetic  energy,  E  is  the  elastic 
energy  and  P  is  the  plastic  energy.  All  energies  are  per 
unit  mass. 

If  fracture  occurs,  a  portion  of  the  elastic  work  is 
converted  into  kinetic  energy  of  the  broken  fragments.  The 
remainder  is  dissipated  at  the  tip  of  the  crack  and  in  running 
the  crack  and  as  surface  energy  or  the  newly  created  fracture 
surfaces.  If  melting  occurs  first,  none  of  the  elastic  work 
is  recovered  as  kinetic  energy,  and  only  plastic  dissipation 
occurs.  Thus,  any  elastic  energy  stored  during  impact  will 
be  recovered  as  kinetic  energy,  converted  to  fracture  surface 
energy  or  dissipated  plastically  or  during  crack  formation 
as  the  penetrator  breaks.  The  energy  balance  after  fracture 
may  be  written 

W  ■  K  +  E*d  (after  fracture)  ,  (63) 

where  a*d  18  the  total  dissipated  energy  per  unit  mass. 

For  many  materials,  EVfd  may  be  thought  of  as  consisting  of 

two  terms:  E,Vdp  which  represents  the  plastic  work  absorbed 

per  unit  mass,  and  E,Vd£  which  is  the  fracture  energy  per 

unit  mass  absorbed  in  the  creation  of  fracture  surfaces.  Thus, 
the  final  energy  balance  may  also  be  written  as 
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(64) 


W  =  K  +  E*dp  +  E*df  . 


In  this  form,  E*dp  includes  the  plastic  dissipation  during 

loading  prior  to  crack  formation  as  well  as  the  energy 
dissipated  to  form  the  crack. 

That  E*d  should  be  an  important  factor  in  penetrator 

behavior  during  impact  is  clear  from  Eq.  (63).  The  work  done  on 
the  penetrator,  if  it  is  not  absorbed  in  E*d  ,  goes  directly 

into  kinetic  energy  of  the  penetrator  in  the  lateral  direction. 
This  leads  to  rapid  spreading  of  the  front  face.  Conversely,  a 
material  with  a  high  E*d  will  absorb  the  work  done  by  the 

the  target  with  less  lateral  spreading  of  the  penetrator. 

To  simulate  the  partitioning  of  energy  in  the  laboratory, 
an  experiment  was  designed  in  which  a  material  specimen  is 
compressed  between  two  anvils  until  fracture  occurs .  Measure¬ 
ment  of  the  time  history  of  the  applied  load  and  the  deforma¬ 
tion  of  the  specimen  yields  the  total  work  done  on  the 
specimen  prior  to  fracture.  Measurement  of  the  kinetic  energy 
associated  with  the  particles  in  the  debris  cloud  following 
comminution  of  the  specimen  yields  a  value  for  K  in  Eq.  (63) . 
The  value  of  E* .  can  then  be  obtained  by  subtracting  K  from 
W  Estimates  of  the  energy  dissipated  during  loading  and 
of  the  fracture  surface  energy  suggest  that  these  two  contri¬ 
butions  to  E,Vj  are  small  and  that  the  primary  energy  dissipa¬ 
tion  mechanism°f or  brittle  materials  is  associated  with  crack 
formation. 


B.  Experimental  Program 

1.  Apparatus  -  Figure  29  is  a  sketch  of  the  EVfd  apparatus. 

The  apparatus  consists  of  a  rigid  frame,  load  application  hardware, 
fragment  trap  hardware,  and  instrumentation. 

The  frame  consists  of  standard  structural  steel  components. 

It  measures  approximately  55  inches  high  by  35  inches  wide  by 
18  inches  deep.  The  load  application  hardware  and  fragment 
trap  hardware  are  mounted  on  and  supported  by  various  platforms 
which  are  fastened  to  the  frame.  To  permit  easy  access  to  the 
test  specimen,  the  frame  rests  on  concrete  blocks. 

A  schematic  illustration  of  the  load  application  hardware 
is  shown  in  Figure  30.  The  material  specimen,  which  in  these 
tests  is  a  spherical  ball,  is  placed  on  the  centerline  of  the 
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apparatus  between  two  crusher  blocks.  The  blocks  are  1-inch 
steel  cubes  which  have  a  Rockwell  hardness  of  C-55.  Since 
compression  of  the  ball  produces  a  small  permanent  indenta¬ 
tion  in  the  blocks,  each  surface  can  be  used  for  only  one 
test. 


Load  is  applied  to  the  ball  using  a  Hanna  MT2  Hydraulic 
Cylinder.  The  cylinder  has  a  bore  diameter  of  4.00  inches  and 
a  rod  diameter  of  1.75  inches.  It  can  generate  a  maximum 
compressive  load  of  approximately  25,000  Ibf.  The  magnitude 
and  rate  of  application  of  the  load  can  be  controlled  by 
regulating  the  hydraulic  pressure  in  the  cylinder.  The 
operating  range  for  the  cylinder  is  200  psig  to  2,000  psig. 

At  the  maximum  cylinder  pressure,  the  rod  velocity  is  approxi¬ 
mately  0.2  inches  /second. 

The  cylinder  is  attached  to  the  upper  beams  of  the  steel 
frame  with  a  trunion  mount.  A  cylinder  stabilizer  platform, 
located  at  the  bottom  of  the  cylinder,  prevents  rotation  of 
the  cylinder  and  maintains  vertical  alignment.  Another 
platform,  the  rod  stabilizer  platform,  is  used  to  maintain 
vertical  alignment  of  the  rod,  A  nylon  phenolic  bushing  in 
the  center  of  the  rod  stabilizer  platform  permits  freedom 
of  movement  for  the  rod.  A  compressive  load  is  applied  to 
the  test  specimen  when  the  rod  moves  downward. 

Load  is  transmitted  from  the  rod  to  the  top  crusher  block 
through  a  steel  rod  cap.  This  cap  has  a  threaded  hole  to 
accept  the  end  of  the  rod.  The  face  of  the  rod  rests  flat  on 
the  inside  surface  of  the  cap.  The  top  crusher  block  is 
ground  flat  and  is  pinned  to  the  rod  cap. 

The  bottom  crusher  block  is  pinned  to  a  steel  anvil 
which  carries  the  load  to  the  load  cell.  The  anvil  has  a 
threaded  male  connection  which  fits  the  load  cell.  The  load 
cell  is,  in  turn,  mounted  on  a  steel  platform  which  is 
fastened  to  the  bottom  beams  of  the  frame. 

Attached  to  the  mid-section  of  the  rod  is  a  collar  and 
fixture  mount  for  the  linear  variable  differential  transformer 
(LVDT)  which  is  used  to  measure  the  displacement  of  the  rod. 

The  LVDT  measures  the  relative  displacement  between  the  rod 
and  the  load  cell  platform. 

An  instrument  plate  is  mounted  to  the  under  side  of  the 
rod  stabilizer  platform.  Four  ballistic  pendulums  are  hung 
from  this  plate.  A  drawing  of  one  of  the  pendulums  1b  shown 
in  Figure  31(a).  Each  pendulum  consists  of  a  1  in^  balsa  wood 
box  with  one  side  removed.  A  styrofoam  block  is  placed  in 
the  box  and  the  front  surface  -  the  surface  exposed  to  the  test 
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specimen  -  is  coated  with  Dow  silicone  grease.  The  mass  of 
the  box  and  styrofoam  is  approximately  2.5  grams.  Prior  to 
impact,  the  front  face  of  the  styrofoam  block  is  adjacent  to 
but  does  not  touch  the  crusher  blocks  (as  shown  in  Figure  30). 
The  four  pendulums  completely  encircle  the  test  specimen  as 
shown  by  the  schematic  in  Figure  31(b).  For  most  of  the  tests, 
the  pendulums  captured  more  than  75%  of  the  mass  of  the  test 
specimen. 

The  arm  of  each  pendulum  is  approximately  four  inches 
long  and  is  made  of  plastic.  For  rigidity,  tne  cross-section 
of  the  arm  is  made  by  gluing  two  I-beam  cross-sections  together 
at  the  flanges.  Each  pendulum  hangs  from  a  separate  shaft 
which  rotates  when  the  pendulum  is  struck.  Each  end  of  the 
shaft  is  mounted  in  a  frictionlesB  pivot  (Bendix  Model  5008-400 
Flexural  Pivot)  which  acts  as  a  torsional  spring.  The  tor¬ 
sional  spring  rate  of  each  pivot  is  6.54  in- lbf /radian.  The 
springs  are  designed  to  produce  a  15  degree  rotation  when 
the  pendulum  is  hit  by  a  1  gram  fragment  traveling  at  75 
feet  per  second.  Also  attached  to  each  shaft  is  a  rotary 
variable  differential  transformer  (RVDT)  to  measure  the  angu¬ 
lar  rotation  of  the  shaft. 

2.  Instrumentation  -  The  instrumentation  consists  of 
one  load  cell .one  lVdt ,  four  RVDT's,  and  one  recording 
device . 

The  load  cell  is  a  Strainsert  Universal  Flat  Load  Cell 
(Model  No.  FL25U-3SG) .  It  is  a  single  bridge  design  with  a 
maximum  capacity  of  25,000  lbf.  The  cell  is  excited  by  a 
12  V  battery  and  produces  a  DC-voltage  output.  It  has  been 
calibrated  against  a  Bourdon  pressure  gauge  which  measures 
the  pressure  in  the  hydraulic  cylinder  and  is  linear  over 
the  full  range  of  application  within  0.5%  of  full  scale. 

The  displacement  of  the  test  specimen  is  measured  with 
a  Trans-Tek  Model  241-000  Displacement  Transducer.  The  trans¬ 
ducer  is  an  integrated  package  consisting  of  a  precision 
linear  variable  differential  transformer,  a  solid  state 
oscillator,  and  a  phase-sensitive  demodulator.  The  LVDT  is  a 
transducer  that  converts  mechanical  displacement  into 
electrical  output.  The  core,  when  displaced  axially  within 
the  coil  assembly,  produces  a  voltage  change  in  the  output 
directly  proportional  to  the  displacement  of  the  core.  This 
model  which  is  excited  by  a  12  V  battery  has  a  working  range 
of  0.1  inches  and  is  accurate  to  0.5%  of  full  scale. 

The  rotation  of  the  shaft  is  measured  by  a  Pickering  Model 
No.  21300  Precision  Rotary  Variable  Differential  Transformer. 
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This  transducer  converts  a  mechanical  angular  displacement 
into  an  electrical  output  by  means  of  an  electrical  input 
carrier.  It  consists  of  a  rotor  assembly  to  which  the  shaft 
is  attached  and  a  stator  assembly  in  which  all  windings  are 
contained.  The  electrical  input  is  a  10  KHz,  3  V  RMS  AC 
signal  which  is  generated  by  a  Hewlett-Packard  Oscillator, 

The  output  signal  is  also  AC  and  is  converted  to  DC  external 
to  the  transducer  using  a  separate  demodulation  circuit. 

The  RVDT  is  linear  over  a  displacement  range  of  i  10*  . 

The  voltage  output  from  each  of  the  six  instruments  is 
applied  to  a  separate  mirror  galvanometer  circuit.  Each 
RVDT  uses  a  Honeywell  M40-120A  galvanometer.  The  load  cell 
uses  an  M40-350A  and  the  LVDT  used  an  M400-120.  The  current 
which  passes  through  the  galvanometer  circuit  is  controlled 
by  external  series  and  shunt  resistances.  The  galvanometers 
are  extremely  sensitive  to  small  changes  in  the  signal  current. 

The  galvanometers  are  contained  in  a  Honeywell  (Model 
906C)  Direct  Recording  Visicorder  Oscillograph.  This  device 
transforms  the  input  signal  into  a  moving  beam  of  light 
through  use  of  the  mirror  galvanometers.  It  uses  an  ultra¬ 
violet  light  source  which  is  focused  by  an  optical  system 
on  recording  paper  which  is  highly  sensitive  to  ultra-violet 
light.  Thus,  data  readout  is  immediately  available.  The 
paper  speed  is  controllable  and,  since  the  device  has  14 
channel  capability,  each  output  trace  can  be  placed  on  the 
same  strip  of  paper. 

3,  Test  procedure  -  Prior  to  a  test,  each  test  &  admen 
is  carefully  weighed  and  measured  and  is  then  placed  at  t-he 
center  of  the  bottom  crusher  block.  A  thin  coat  of  Dovr 
silicone  grease  is  applied  to  one  surface  of  each  styrofoam 
block,  ana  the  block  is  then  weighed.  The  rod  Is  positioned 
such  that  the  top  crusher  block  is  approximately  one  inch 
above  the  test  specimen  and  the  pressure  in  the  cylinder  is 
set  to  give  the  desired  loading  rate. 

The  oscillograph  is  actuated  at  the  instant  the  cylinder 
rod  begins  moving.  Note  that  the  top  crusher  block  is  moving 
downward  when  it  comes  in  contact  with  the  test  specimen. 

Load  is  applied  until  the  ball  fractures  or  the  rod  bottans- 
out . 


After  the  test,  the  styrofoam  blocks  are  carefully 
removed  and  weighed  to  determine  the  mass  of  the  pieces  which 
impacted  the  blocks.  Any  debris  which  remains  on  the  crusher 
blocks  is  also  collected  and  weighed.  Finally,  the  crusher 
blocks  are  removed,  and  the  depth  of  the  indentation  is 
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measured  using  a  dial  indicator  gauge. 


4.  Data  analysis  -  The  oscillograph  provides  a  continu¬ 
ous  time  history  of  the  load  and  deflection.  Recall  that 
the  LVDT  shown  in  Figure  30  actually  measures  the  position  of  the 

bottom  of  the  rod  relative  to  the  top  of  the  load  cell  plat¬ 
form.  Therefore,  the  deflection  measurement  includes  the 
indentation  of  the  crusher  blocks.  (The  displacement  of  the 
rod  cap,  anvil  and  load  cell  are  negligible.)  The  area 
beneath  the  force-deflection  curve  up  to  rupture  represents 
the  total  work  done  on  the  ball  and  the  crusher  blocks. 

Equation  (62)  which  describes  the  energy  balance  based  on 
a  rigid  crusher  block  must  be  modified  to  account  for  energy 
absorption  by  the  blocks.  Both  elastic  and  plastic  work  is 
done  on  the  blocks.  It  is  assumed  that  the  stored  elastic 
energy  is  returned  to  the  ball  when  the  ball  fractures. 

However,  the  plastic  work  is  dissipated  in  the  blocks.  Hence, 
for  a  non-rigid  crusher  block  the  energy  equation  may  be 
written 

W  -  K  +  E*d  +  Edb  ,  (65) 

where  E..  represents  the  energy  dissipation  in  the  crusher 
blocks. 

The  quantity  Edb  can  be  estimated  by  testing  a  rigid 

bdll,  Figure  32  shows  the  load-deflection  response  of  the 
system  when  a  tungsten  carbide  ball  1b  compressed  through 
five  cycles  of  loading  and  unloading.  At  the  end  of  each 
cycle  there  is  no  measureable  plastic  deformation  of  the  ball 
so  that  the  measured  def  lectif  on  represents  the  total  indenta¬ 
tion  of  the  crusher  blocks.  For  each  cycle,  the  area  beneath 
the  loading  curve  is  the  total  work  done  on  the  system;  the 
area  beneath  the  unloading  curve  is  the  recovered  elastic 
energy.  The  difference  between  the  two  represents  the  energy 
dissipation  in  the  blocks.  Figure  33  shows  this  result  as  a 
function  of  the  depth  of  the  indentation. 

The  oscillograph  trace  of  the  RVDT  response  is  essentially 
a  step  pulse  followed  by  a  damped  oscillation.  The  magnitude 
of  the  pulse  which  is  calibrated  to  the  rotation  of  the  pen¬ 
dulum  can  be  used  to  estimate  the  kinetic  energy  of  the 
impacting  particles,  It  is  assumed  that  all  of  the  pieces 
which  impact  a  given  pendulum  can  be  lumped  together  to  form 
an  average  particle  which  has  a  mass,  mproj ,  equal  to  the 

total  mass  of  all  the  pieces  and  a  translational  velocity, 

Vproj  ’  norinal  t0  fflce  of  the  pendulum.  The  lumping  of  all 
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of  the  debris  into  one  average  piece  is  admittedly  a  simplifi¬ 
cation.  However,  the  error  is  minimized  by  removing  from 
consideration  those  tests  in  which  the  RVDT  loading  pulse  was 
jagged.  Neglect  of  the  tangential  component  of  the  impact 
velocity  introduces  an  error  of  approximately  15%  in  the 
estimate  for  Vproj  •  This  error  was  estimated  based  on  an 

angle  of  approximately  45°  between  the  center  of  the  ball  and 
the  side  of  the  styrofoam  block  and  an  angle  of  less  than 
15°  to  the  top  of  the  block, 


The  kinetic  energy  of  the  debris  can  be  obtained  from 
the  conservation  equations  as  follows.  The  momentum  equation 
may  be  written 


mproj  Vproj  *  ^mproj+mbox^  Vbox  ’ 

where  m^ox  is  the  mass  of  the  pendulum  and  is  the 

initial  velocity  of  the  pendulun.  The  energy  equation  may 
be  written 
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where  the  first  term  on  the  right  accounts  for  the  change  in 
potential  energy  as  the  pendulum  swings  through  an  arc  of  6 
degrees  (L  is  the  length  of  the  pendulum),  The  second  term 
represents  the  torsional  energy  stored  in  the  springs  of  the 
flexural  pivots.  Frictional  losses  in  the  system  are  accounted 
for  by  using  a  calibration  value  for  kg  .  This  value  is 

within  4%  of  the  nominal  spring  constant  for  each  pendulum. 

For  small  rotations  (i.e.,  less  than  10°)  the  potential 
energy  term  in  Eq.  (67)  can  be  ignored.  Thus, 
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The  initial  velocity  of  the  pendulum  can  be  obtained  from  the 
measured  mass,  the  angular  displacement  of  the  pendulum  and 
the  Known  spring  constant.  Thus  ^proj  can  be  obtained 

from  Lq .  (66)  and  K  ,  the  kinetic  energy  of  the  particles, can 
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be  obtained  from 


proj  proj 


5,  Test  program  -  Table  4  summarizes  the  materials  which 
were  investigated  and  the  pertinent  results.  A  detailed  dis¬ 
cussion  of  the  results  is  presented  in  the  next  section. 

C.  Test  Results 

1.  Borosilicate  (Pyrex)  glass  -  Figures  34  -  36  show  the 
force-deflection  data  for  borosilicate  glass  (Corning  7740 
polished  pyrex)  spheres  for  three  strain  rate  regimes.  For 
the  purposes  of  tnis  report,  the  strain  rate  is  defined  as 
an  average  value  for  the  entire  test 


where 


is  the  measured  displacement, 


is  the  initial 


diameter  of  the  ball  and 


is  the  time  to  rupture.  Figure 


34  shows  the  data  for  four  tests  which  were  performed  at 
approximately  the  same  strain  rate.  This  figure  suggests 
the  degree  of  repeatability  of  the  data.  Three  of  the  four 
loading  curves  virtually  overlap.  However,  the  rupture  load 
differs  by  approximately  350  lbf ,  and  the  deflection  at 
rupture  by  approximately  .002  inches  between  the  tests.  For 
reference  purposes,  the  classical  solution  to  the  problem 
is  also  shown.  This  solution  which  was  first  derived  by 
Hertz  can  be  found  in  Ref.  7. 

Figure  35  shows  the  loading  data  for  tests  at  three 
intermediate  strain  rates,  and  Figure  36  shows  the  data  for 
the  highest  strain  rates  which  were  achieved.  The  rupture 
load  and  the  deflection  at  rupture  for  each  of  these  tests  is 
shown  in  Figure  37.  Although  there  is  considerable  scatter 


7.  Bergstrom,  B.  H.  and  Sollenberger ,  C.  L. ;  Kinetic  Energy 
Effect  in  Single  Particle  Crushing!  Presented  at  the  AIME 
Meeting,  St.  Louis,  February  1961. 
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in  the  data,  there  does  appear  to  be  a  discernable  trend 
toward  increasing  rupture  load  and  deflection  at  the  higher 
strain  rates. 

The  area  beneath  the  loading  curve  represents  the  total 
work  done  on  the  test  specimen  and  crusher  blocks.  Figure 
38(a)  shows  these  results  and,  again,  there  appears  to  be  a 
tendency  for  W  to  increase  with  strain  rate.  Where  does 
this  energy  go  when  the  ball  is  fractured?  Figure  38(b) shows 
the  measured  kinetic  energy  and  Figure  38(c)  shows  the  ratio 
of  the  kinetic  energy  to  the  work  done.  For  this  material 
and  these  strain  rates,  approximately  15%  of  the  work  is 
recovered  as  non-dissipative  kinetic  energy.  It  is  important 
to  note  that  this  percentage  appears  to  be  roughly  constant 
for  this  range  of  strain  rates.  Keep  in  mind  that  this 
percentage  is  probably  somewhat  low  for  the  reasons  discussed 
earlier,  and  also  because  the  kinetic  energy  associated  with 
particles  which  rebound  between  the  crusher  blocks  is  not 
measured.  The  mass  not  captured  by  the  pendulums  is  generally 
less  than  207.  of  the  ball  mass  so  this  error  cannot  be  much 
larger  than  20%.  Even  if  the  kinetic  energy  measurement  is 
in  error  by  a  factor  of  two,  however,  it  is  clear  that  more 
than  70%  of  the  work  is  dissipated. 

Energy  can  be  dissipated  in  various  modes  such  as  heat 
or  in  free  surface  energy  or  to  initiate  and  run  cracks. 
Estimates  for  the  energy  dissipation  via  some  of  these 
modes  can  be  made.  Figure  39  shows  one  cycle  for  a  test 
in  which  the  load  was  released  prior  to  fracture.  The 
unloading  curve  overlaps  the  loading  curve.  Hence,  there  is 
no  plastic  dissipation  in  the  ball  prior  to  fracture.  The 
indentation  in  the  crusher  blocks  was  less  than  .001  inches. 
Hence,  the  dissipation  in  the  crusher  blocks  was  less  than 
10%  of  the  total  work.  An  estimate  was  also  made  of  the 
energy  associated  with  the  surface  area  created  during  the 
fracture.  A  test  was  conducted  in  which  all  of  the  broken 
pieces  were  collected,  and  an  estimate  was  made  of  the  total 
surface  area.  Based  on  a  theoretical  value  for  the  free 
surface  energy  per  unit  area  for  glass  of  1.7  x  103  ergs/cm^, 
it  was  estimated  that  less  than  1%  of  the  total  work  was 
converted  to  free  surface  energy.  Thus,  nearly  all  of  the 
dissipated  energy  was  associated  with  the  formation  and 
running  of  the  cracks, 

2.  Soda  lime  glass  -  The  loading  curves  for  soda  lime 
glass  (dQ  -  .25  inches)  are  shown  in  Figures  40  and  41  for 

two  strain  rate  regimes.  The  rupture  load  and  the  total 
deflection  at  rupture  for  these  tests  is  summarized  in 
Figure  42.  In  general,  the  rupture  load  iE  considerably 
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higher  for  this  glass  than  for  the  borosilicate  glass, 
although  the  deflections  at  rupture  are  comparable ,  The 
partitioning  of  energy  is  shown  in  Figure  43.  It  is  clear 

that  less  than  15%  of  the  total  work  was  recovered  as 
non-dissipative  energy,  Thus,  nearly  85%  of  the  energy  was 
dissipated.  Again,  the  plastic  dissipation  during  deforma¬ 
tion,  the  plastic  dissipation  in  the  crusher  blocks,  and  the 
fracture  surface  energy  were  a  small  percentage  of  the  total 
energy  dissipation.  The  bulk  of  the  energy  was  dissipated 
during  the  formation  and  running  of  the  cracks.  Also  note 
that  the  ratio  of  K/W  is  approximately  constant  for  this 
range  of  strain  rates. 

A  series  of  tests  was  also  conducted  on  larger  soda 
lime  glass  balls  (d  -  0.5  inches).  The  loading  data  for 
these  tests  are  shown  in  Figure  44.  There  appears  to  be  a 
definite  tendency  for  the  rupture  load  to  increase  with 
strain  rate.  Note  that  these  loads  are  roughly  four  times 
the  rupture  loads  for  the  0.25  inch  balls  so  that  the 
rupture  load  approximately  scales  with  the  area.  In  addition, 
the  displacements  are  roughly  twice  the  values  for  the  0,25 
inch  balls  so  that  displacement  roughly  scales  with  diameter. 

In  these  tests,  the  non-dissipated  kinetic  energy  which 
was  recovered  varied  between  22%  and  387,  of  the  total  stored 
energy.  Thus,  the  dissipated  energy  varied  between  roughly 
60%  and  80%  of  the  stored  energy.  Again,  this  dissipated 
energy  could  only  be  associated  with  the  formation  and 
running  of  the  cracks. 

3.  Aluminum  oxide  -  Five  tests  were  conducted  on 
alumina  (99 . 9%  pure')  balls .  The  loading  data  are  shown  in 
Figure  45,  There  was  approximately  a  20*4  increase  in  the 
rupture  load  and  the  displacement  at  rupture  for  this  range 
of  strain  rates.  The  recovered  kinetic  energy  never  exceeded 
4%  of  the  total  stored  energy.  Thus,  the  dissipated  energy 
was  at  least  96%  of  the  stored  elastic  energy. 

4.  Tungsten  carbide  -  The  tests  on  tungsten  carbide 
balls  also -showed  a  very  small  nondissipative  energy  component. 
The  loading  curve  for  one  of  the  tests  is  shown  in  Figure  46, 
The  work  done  on  the  ball  during  this  test  was  approximately 
112  in-lbf.  The  total  indentation  in  the  crusher  blocks  was 
approximately  0,019  inches  and,  using  Figure  33,  approximately 
80  in-lbf  were  dissipated  in  the  blocks.  Approximately  50% 

of  the  ball  was  captured  in  the  pendulums,  and  the  kinetic 
energy  of  these  particles  was  less  than  1  in-lbf.  Thus,  of 
the  32  in-lbf  of  work  done  on  the  ball,  more  than  97%  was 
dissipated  during  fracture. 
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5.  Metals  -  A  number  of  tests  were  conducted  on  metal 
balls'!  For  the  soft  metals  such  as  1013  steel,  1100-F 
aluminum,  and  2717-T4  aluminum  the  ball  failure  was  ductile, 
i.e. ,  the  ball  was  compressed  into  a  flat  disk  and  did  not 
fracture,  Hence,  all  of  the  energy  was  plastically  dissipated 
during  deformation, 

However,  the  hard  chrome  steel  balls  (AISI  52100)  did 
exhibit  brittle  fracture.  Figure  47  shows  the  loading  data 
for  one  of  the  tests  conducted  with  this  material.  The 
rupture  load  was  approximately  13,000  lbf.  The  total  work 
was  approximately  3/0  in-lbf.  Of  this  total,  approximately 
1  in-lof  was  recovered  as  nondlssipated  kinetic  energy. 

Thus,  virtually  all  of  the  work  done  on  the  ball  was  dissipated. 

6.  Plastics  -  Tests  were  conducted  on  several  plastic 
balls  including  polycarbonate  (Lexan) ,  acrylic  (Lucite) , 
and  phenolic.  For  the  polycarbonate,  the  loading  rates 
were  too  low  to  produce  brittle  fracture;  the  ball  was  com¬ 
pressed  into  a  flat  disk. 

Figure  48  shows  the  loading  data  for  the  phenolic  ball. 

For  this  test,  the  total  work  was  approximately  5  in-lbf.  The 
recovered  kinetic  energy  was  only  0.13  in-lbf.  Thus  more  than 
97%  of  the  energy  was  dissipated. 

Figure  49  shows  the  loading  data  for  the  acrylic  ball. 

The  results  for  this  material  were  similar  to  those  for  the 
phenolic.  The  total  work  was  approximately  4.3  in-lbf  and 
less  than  1%  of  this  energy  was  recovered  as  kinetic  energy. 

The  remainder  was  dissipated. 

A  summary  of  the  results  of  these  tests  is  given  in 
Table  4.  For  the  range  of  strain  rates  tested,  the  dissipative 
energy  E,V(^  is  always  much  larger  than  the  recoverable  kinetic 

energy.  The  primary  component  of  the  dissipated  energy  appears 
to  be  the  energy  associated  with  the  formation  and  running  of 
cracks  when  the  ball  fractures.  Extrapolation  of  these  results 
to  impact  strain  rates  is  difficult,  and  it  is  not  clear  that 
the  partitioning  obtained  in  the  laboratory  will  occur  during 
impact.  However,  it  is  clear  that  an  efficient  penetrator 
materia,  must  be  able  to  dissipate  most  of  the  work  done  on 
it  by  the  target;  the  recoverable  kinetic  energy  must  be  kept 
to  a  minimum. 
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VI .  CONCLUSIONS 


1.  A  wide  range  of  target  materials  has  been  tested  in  the 
A.R.A.P.  Impact  Facility,  and  their  impact  properties  have 
been  evaluated  using  the  integral  theory  of  impact.  The 
materials  tested  include  soft  and  hard  metals,  ceramics, 
metal  alloys,  plastics  and  composites.  The  results  demon¬ 
strate  that  the  impact  performance  of  a  target  material  can 
be  characterized  using  the  two  material  properties  E*„  and 

E*e  •  P 

2.  A  theory  has  been  developed  which  relates  and  E*g 

to  more  fundamental  material  properties.  The  theory  yields 
simple  equations  which  can  be  used  as  a  preliminary  screening 
tool  to  suggest  candidate  materials  for  particular  Impact 
applications. 

3.  A  long  rod  penetrator  model  has  been  developed.  The 
model  contains  two  cells,  satisfies  global  conservation  equa¬ 
tions,  and  accounts  for  material  strength  using  the  A.R.A.P. 
concept  of  adiabatic  hardness  pE*  .  The  numerical  code 
based  on  this  model  predicts  penetration  depth,  ballistic 
limit  velocity,  residual  mass  and  residual  velocity  in  good 
agreement  with  data  for  the  high  velocity  impact  of  homogeneous 
rods. 

4.  A  fundamental  experiment  has  been  performed  at  laboratory 
strain  rates  which  demonstrates  that  most  of  the  strain  energy 
stored  in  a  penetrator  during  compression  is  dissipated 

when  the  penetrator  fractures.  The  primary  dissipation 
mechanism  appears  to  be  associated  with  the  formation  and 
running  of  cracks. 
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APPENDIX  A 


MATERIAL  QUALIFICATION  TEST  DATA 

This  appendix  contains  a  summary  of  the  data  obtained 
during  the  material  qualification  tests.  For  each  material, 
the  following  items  are  tabulated: 


(1) 

Test  number 

(2) 

Target  identification  number 

(3) 

Target  thickness,  inches 

(4) 

Projectile  material,  WC  • 

ST  - 
STB  - 

tungsten  carbide 
chrome  steel 
soft  steel 

(5) 

Projectile  mass,  grams 

(6) 

Projectile  diameter,  inches 

(7) 

Projectile  velocity,  feet /second 

(8) 

Crater  depth,  inches 

(9) 

Note  for  test  anomalies 

A  -  Projectile  remained  embedded  in  target 
B  ■  Projectile  broken 
C  ■  Target  backface  cracked 

D  =  Penetration  greater  than  half  the  target 
thicknesB 

E  ■  Target  backed  with  rubber. 

Note  that  the  data  are  tabulated  in  both  metric  and  English 
standard  units.  The  conversion  factors  between  the  systems 
are 

1  gram  «*  2.205  x  10”^  lbm 
1  meter  *  3 . 281  feet 
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APPENDIX  B 


E*  EVALUATION 

The  results  of  the  E*  evaluations  for  each  of  the 
target  materials  which  were  qualified  are  shown  in  Figures 
B-l  to  B-14.  The  procedure  by  which  these  values  were  obtained 
was  described  in  Chapter  II.  Briefly,  these  values  of  E* 
when  used  in  the  Integral  theory  result  in  a  computed  penetra¬ 
tion  which  matches  the  measured  penetration  for  a  given  impact 
velocity.  The  data  are  plotted  versus  nondimensional  penetra¬ 
tion  depth  because  of  the  functional  dependence  of  E*  on 
p/d  which  was  noted  in  Chapters  II  and  III. 
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Figure  B-9 
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APPENDIX  C 

USER'S  GUIDE  TO  ROD  CODE 

This  appendix  is  a  user's  guide  for  the  ROD  code.  It 

contains : 

(1)  a  summary  of  ths  system  of  nondimens tonal  equations; 

(2)  an  explanation  of  the  method  of  solution; 

(3)  FORTRAN  listings  of  the  mainline  program  and  all 
subroutines ; 

(A)  input  specifications; 

(5)  output  translation;  and 

(6)  a  sample  case. 

A.  Nondlmensional  Equations 

The  system  of  equations  which  defines  the  penetration  of 
a  long  rod  projectile  is  summarized  in  Table  3.  This  system 
consists  of  14  first-order  ordinary  differential  equations. 

It  is  convenient  to  solve  these  equations  in  nondlmensional 
form.  The  nondlmensional  parameters  are  summerlased  in 
Table  C-l. 

The  solution  technique  requires  that  the  differential 
equations  be  in  the  form 

V  -  g  -  F(x,Y)  ,  (C.l) 

where  Y  denotes  a  column  vector  of  the  dependent  variables, 
and  x  la  the  independent  (time  variable).  Without  going 
into  the  details  of  the  algebraic  manipulations,  it  is 
possible  to  rearrange  the  equations  to  obtain  the  system 
summarized  in  Table  C-2.  In  the  order  in  which  the  equations 
are  written,  each  derivative  is  expressed  in  terms  of  the 
Independent  variable  and  the  preceding  dependent  variables 
as  required  by  Eq.  (C.l). 


TABLE  C-l 


NORMALIZATION  PARAMETERS 


Quantity  Dimensional  Nondimens Iona 1 

Symbol  Symbol 


Shaft  Velocity 

V. 

m 

v,/v+ 

Head  Center  of  Mass  Velocity 

Vcm 

Ycm" 

v„m/v+ 
cm  _l 

Head  Front  Face  Velocity 

Vi 

Yi 

m 

yv; 

Head  Edge  Velocity 

VT 

Yt 

u 

Vvt 

Interface  Velocity 

\ 

vt 

m 

Vb/V+ 

Mass  in  shaft 

ms 

m 

“s 

m 

ra8/mp 

Mass  in  head 

m 

m 

"h/mP 

Mass  eroded 

®b 

% 

m 

Radius  of  head 

b 

m 

b 

m 

h/,lv 

Half  thickness  of  head 

H 

i 

m 

ih 

Depth  of  Penetration 

P 

P 

m 

P/"P 

Total  Kinetic  Energy 

K 

K 

m 

K/K 

Shaft  Kinetic  Energy 

Ks 

m 

kA 

K  /K 
cm  p 

Head  Center  of  Mass  Kinetic  Energy 

^cm 

K  * 

cm 

Head  Relative  Kinetic  Energy 

Kr 

K 

r 

M 

K  /K 
r  p 

Density  Ratio 

m 

P 

m 

Pt/c 

Rod  L/D 

Hydrodynamic  Mode  Energy  of  Target 

Rod  Energy  Dissipation 

Adiabatic  Yield  Strength  of  Rod 

E* 

E*d 

0 

L  ■ 

l*dm 

3  ■ 

E*/V+2 
oh  PV+ 

Time 

t 

T 

- 

v+  , 

IT  t 

Normalization  Parameters  are 

P 

V+  -  characteristic  velocity 

O 

o 

o 

cH 

■c 

ft/sec 

m  »  initial  mass  of  rod  -  p 

2p  *  initial  length  of  rod 

l  IT 

Kp  -  characteristic  energy  - 

1_  V4-2 

7mpv 

a  -  radius  of  rod 
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TABLE  C-2 

NONDIMENS I ONAL  SYSTEM  OF  EQUATIONS 


V. 


®b 

% 

^cm 


K 

^cm 


k; 


K 


•  t 

V, 


V 

VT 


a  f 

r'  . 3 

16  L2  b  2  Vt  fh 


6  f 


- 1 


•  “b 

u. 


±  r 
*.  L 


^cm®b *.  K  -  Vcmm  i  -  4b2L2'p(i*  +  ^Vi2) 


v/  «;  +  2  a,  v,  v' 

5cm^  +  2Sh5Cm  vi« 


+  2  o  £ 


(v  -  V  ^ 

\  8  cny 


«i  [v? + H’*  •  Wl  £h + 8Pe  1 L  ( %  -  vi)  (** +  -t  O 

m.  L  ~r~  -E* 


Kl  +  Km  +  K 

s  cm  r 

(2  ■  fh) 

51  -  5cm 

V.  I  2*B  •  3M 

.CTn._  4fi*  - 

C  -  EE  c 

2B 


wher*  f  ■  0  mg  ■  0 
£g  -  !  mg  >  0 


fh"0  b<<1  +  eo>bo 
fh"l  6i<1  +  eo)bo 


L-  2p/2a 

«  -  VB  -  Vb 

-  «  v,  -  V 
E  Vl  cm 


A'r13Vfi4 

L 

D-l/b 

E  -  1/2  -  £  /Sh 


;*t2  +  (vi-v0 J 


C-( 
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cm 


mh 


b'V^  l 


^+T(b'+vt) 


aa<tir>ff‘‘tffeSaAyjaytfiia^aL  r,  takua&ui «*:,  vu'juliuk  * 


B.  Method  of  Solution 


The  numerical  problem  is  to  solve  a  system  of  simultane¬ 
ous  first-order  ordinary  differential  equations  given  the 
initial  value  for  each  variable.  The  solution  technique  uses 
Hamming's  modified  predictor-corrector  method  for  the  solution 
of  general  initial  value  problems.  This  method  is  a  stable, 
fourth-order  integration  procedure  that  requires  the  evalua¬ 
tion  of  the  right-hand  side  of  the  system  only  two  times  per 
step.  This  is  a  great  advantage  over  other  methods  with  the 
same  order  of  accuracy,  such  as  the  Runge-Kutta  method,  which 
requires  four  evaluations  of  the  right-hand  side  per  step. 
However,  the  method  is  not  self- starting;  that  is,  the 
functional  values  at  a  single  previous  point  are  not  enough 
to  get  the  functional  values  ahead.  Therefore,  to  get  the 
starting  values,  a  special  Runge-Kutta  procedure  followed  by 
one  iteration  step  is  added  to  the  predictor- corrector  method. 
In  the  ROD  code,  the  subroutine  which  performs  the  integration 
is  HFCG  -  a  listing  of  which  is  provided  later  in  this 
appendix.  The  routine  is  part  of  IBM's  scientific  software 
package,  and  precise  details  of  the  numerical  procedure  are 
available  (Ref.  8). 

The  ROD  Code  is  run  in  discrete  time  steps;  the  size  of 
the  steps  is  controlled  by  an  input  parameter.  The  predictor- 
corrector  method  is  used  to  solve  an  initial  value  problem  for 
each  time  interval;  the  solution  at  the  end  of  a  given  time 
step  is  the  initial  condition  for  the  next  stap. 

C.  ROD  Code  Listings 

A  FORTRAN  source  listing  of  the  ROD  code  is  provided 
at  the  rear  of  this  appendix.  The  code  consists  of  a  mainline 
routine  RODML  and  four  subroutines  RODDE,  RODES,  HPCG,  and 
OUTP.  These  routines  perform  the  following  functions. 

RODML  -  mainline  routine  which  controls : 

(1)  input  and  output, 

(2)  variable  initialization, 

(3)  integration  of  aquations, 

(4)  computation  of  auxiliary  quantities,  and 

(5)  run  termination. 


Ref.  8>  IBM  large  computer  scientific  subroutine  package. 
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RODES  -  subroutine  which  evaluates 

HPCG  -  subroutine  which  integrates  the  system  of 
equations 

RODDE  -  subroutine  which  evaluates  the  derivatives 
as  required  by  HPCG 

OUTP  -  external  output  subroutine  required  by 

HPCG.  It  is  not  used  in  ROD  except  insofar 
as  to  satisfy  the  EXTERNAL  statement. 

D.  Input  Specifications 

The  ROD  code  can  be  run  using  card  input.  The  input 
specifications  are  shown  below.  FORTRAN  format  statements 
are  shown  in  parentheses. 

CARD  1  (115)  Number  of  Cases 

NCASE  -  number  of  impact  cases  which  will  be  run 
(no  limit) 

For  each  case,  the  following  cards  are  required. 

CARD  2  (40A2)  Title  Card 

ITITL  -  title  or  heading  which  will  appear  on  each 
page  of  print-out 

CARD  3  (6F10.0)  Rod  Specifications 

VELOC  -  normalized  Impact  velocity,  equal  to  V/1000 
where  V  is  the  impact  velocity  in  ft/ sec 

ELOVD  -  length  to  diameter  ratio  of  the  rod 

EFREE  -  normalized  E*j  of  the  rod  material,  equal 
to  .01052  x  E*  where  E*  has  the  units 
Btu/lbm 

SIGMA  -  normalized  adiabatic  yield  stress  of  the 
rod,  equal  to  .01052  x  E*  where  E*  has 
the  unitB  Btu/lbm 

EPSO  -  the  plastic  strain  permitted  in  the  head 
before  erosion  occurs,  present  value  0.36 

THICK  -  normalized  thickness  of  the  target,  equal 
to  T/ip  where  T  is  the  thickness  of 


the  target  and  ft  is  the  length  of  the 
projectile  p 

OBLIQ  -  angle  of  incidence,  in  degrees,  from  the 
normal 

CARD  4  (215)  Flag  Specifications 

NPRNT  -  number  of  integration  steps  between  print¬ 
outs  (usually  10) 

NLAYR  -  number  of  separate  material  layers  in  the 
target  (maximum  value  -  5) 

CARDS  5-9  (4F10.0)  Target  Specifications  , 

A  separate  card  is  needed  for  each  layer  of  the  target. 

-  ratio  of  layer  material  density  to  pro¬ 
jectile  material  density 

-  normalized  hydrodynamic  mode  energy,  E*  , 
of  the  layer,  equal  to  .02505  x  E*  and 
E*  has  the  units  Btu/lbm 

-  the  normalized  line-of- sight  thickness  of 
the  layer,  equal  to  Tp/ftp  where  Tm  is 
the  line-of-slght  thickness  of  the  layer 
and  ftp  is  the  length  of  the  projectile. 
Note  that 

NLAYR 

THICK  -  COS  (OBLIQ)  x  T*  Tm/  ft  . 

S3 

-  the  integration  step  size  in  nondimens ional. 
time,  usually  0.002.  Recall,  the  non- 
dimensional  time,  t  ,  is  relatrd  to  real 
time,  t  ,  by 

t  ■  1000  t/ ft_  . 

P 

The  cards  for  all  impact  cases  must  be  input  prior  to 
execution  of  the  first  case.  The  code  will  process  each  case 
sequentially  provided  the  preceding  case  is  not  terminated  by 
an  error.  If  an  error  occurs,  the  execution  of  all  subsequent 
cases  is  terminated.  A  successful  termination  of  the  case 
occurs  when  one  of  the  following  occurs: 


DENS 

ESTR 

THIK 

DTAUI 
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(1)  the  total  kinetic  energy  of  the  rod  is  less  than 
0.17„  of  the  initial  kinetic  energy; 

(2)  the  axial  momentum  of  the  rod  is  less  than  0.1%  of 

the  initial  axial  momentum; 

(3)  the  front  face  of  the  head  of  the  rod  penetrates  a 
distance  equal  to  the  total  thickness  of  the  target. 

There  are  various  conditions  for  which  abrupt  termination 
of  the  code  will  occur.  When  the  program  terminates  prior  to 
successful  completion  of  a  case,  an  error  message  is  printed 
These  errors  include: 

ERROR  1  -  denotes  a  situation  which  yields  an  imaginary 
solution  of  the  equations.  This  error  is  the  result  of  global 
modeling  of  the  deformation  field  in  the  head.  There  is  no 
solution  to  the  problem  using  this  version  of  the  ROD  code 
although  a  future  version  which  utilizes  a  somewhat  different 
formulation  for  the  system  of  equations  will  eliminate  the 
error. 

ERROR  2  -  denotes  an  error  in  the  solution  subroutine 
HPCG.  The  predictor-corrector  algorithm  cannot  integrate  the 
equations  within  its  own  Internal  tolerances.  The  problem 
may  be  solved  by  changing  the  time  step  -  which  is  an  input 
item  -  or  the  internal  tolerances  within  HPCG  -  which  requires 
code  editing. 

ERROR  3  -  denotes  excessive  build-up  of  integration  error. 
Try  a  smaller  time  step. 

ERROR  4  -  denotes  that  the  code  requires  more  than  2,000 
integration  steps  to  complete  the  case.  This  error  is 
included  to  prevent  looping.  Try  running  with  a  larger  time 
step. 

In  addition  to  these  error  messages,  there  are  two 
additional  messages  which  will  occur  which  are  not  errors. 

END  OF  LAYER  -  denotes  that  the  front  face  of  the  pro¬ 
jectile  has  penetrated  a  distance  equivalent  to  the  thickness 
of  a  layer. 

END  OF  SHAFT  -  denotes  that  the  shaft  has  been  completely 
eroded,  and  only  the  head  of  the  rod  remains. 


E.  Sample  Case  and  Output  Translation 


A  sample  case  is  presented  to  Illustrate  the  input 
specifications  and  to  translate  the  output.  The  case  is  the 
normal  impact  of  a  tungsten  alloy  rod  (Mallory  3000)  into  a 
semi- infinite  target  or  rolled  homogeneous  steel  armor.  The 
impact  velocity  is  5,046  feet  per  second.  The  value  of  E* 
for  the  rod  material  is  50  Btu/lbm  and  for  the  RHA  target 
E*  is  200  Btu/lbm.  This  case  corresponds  to  an  actual  test 
performed  by  BRL  for  ARAP. 

The  input  cards  for  this  case  are  as  follows: 


CARD  1  1, 

CARD  2  MALLORY  3000  ROD  INTO  RHA  AT  5046  FPS 

CARD  3  5.046,10.1, .526, .526, .36,1000. ,0. 

CARD  4  10,1, 

CARD  5  0.462,5.01,1000. , .002, 

An  explanation  for  each  of  these  Inputs  is  given  below. 


CARD  1  NCASE  ■  1  One  Impact  case  will  be  calculated. 

CARD  2  Title  and  page  heading  for  case 


CARD  3  VELOC 

ELOVD 

EFREE 


V  .  5*046  fps  .  5  046 

V+  1 , 000  fps 

£  80.26  mm 

_E  - - -  10.1 

D  7, 92  mm 


E*d  -42E* 
J+i  *  !pr 


2.503x10*1^ 
2^6  ftz  Btu/lbm 


sec 

.01052  x  50  -  526 

Note  the  conversion  factor  required  to 
make  EFREE  dimensionless. 


a  . 42p  E* 

_ _  «  _ ~  P- 

PpV+7  PpV+2 


.42  E* 


SIGMA 


EFREE  -  .526 


EPSO  -  .36  (model  parameter) 

THICK  *  T/£  *  1000  (very  large  value  to 

approximate  semi-infinite  target) 

OBLIQ  -  0.  (normal  impact) 

CARD  4  NPRNT  «  10  The  solution  for  every  tenth  inte¬ 
gration  step  will  be  printed.  In 
addition,  printouts  automatically  occur 
for  the  first  step  and  for  significant 
steps  such  as  layer  penetration. 

NLAYR  -  1  The  target  consists  of  one  layer. 


CARD  5  DENS  -  -=■ 


7.86  gm/cc 


-  .462 


ESTR  - 


Pp  17.0  gm/cc 

E*  200  Srjl-  a  ft2/ sec2 

*  .  x  2.505  xlO4 - -5.01 

v+2  106  £t±  Btu/lbn 


sec 


THIK  -  -SJ.  -  1000  (very  large  value  to  approxi¬ 


mate  semi-infinite  layer) 


DTAUI  -  .002  nondimensional  time  step. 


The  output  for  this  case  is  shown  on  the  following  pages . 
The  top  of  the  first  output  page  contains  a  summary  of  the 
input  specifications  for  the  case.  This  is  followed  by  three 
lineB  of  headings  for  the  output  parameters,  Table  C-3  pro¬ 
vides  a  translation  for  each  of  these  parameters.  Each  three 
line  grouping  beneath  these  headings  provides  the  solution  of 
the  svstem  of  equations  at  the  indicated  time.  Thus,  the 
time  history  (or  the  history  versus  depth  of  penetration)  for 
any  of  the  variables  is  readily  available.  The  final  grouping 
of  three  lines  is  the  soluticn  at  the  final  time  step.  In  this 
case,  the  penetration  terminates  at  t  -  .402  (which  corre¬ 
sponds  to  t  *  106  y  sec), and  the  rod  has  penetrated  a  distance 
PEN  -  .954  (which  corresponds  to  ,954  x  80.26  mm  -  76.6  mm). 

The  run  is  terminated  at  this  point  because  the  total  kinetic 
energy  in  the  head  is  less  than  0.1%  of  the  initial  kinetic 
energy  in  the  rod. 
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Table  C-3 


OUTPUT  TRANSLATION 


Output  Heading 

Symbol 

Definition 

TAU 

T 

Table  C-l 

VS 

Table  C-l 

VCM 

!« 

Table  C-l 

V-PERP 

Vi 

Table  C-l 

VT 

Yt 

Table  C-l 

VB 

Vb 

Table  C-l 

WIDTH 

-  2  x  b  (Table  C-l) 

THICK 

- 

•2x1  (Table  C-l) 

PEN 

•B 

P 

Table  C-l 

MH 

*h 

Table  C-l 

MS 

5. 

Table  C-l 

MB 

°b 

Table  C-l 

MASS  RATIO 

■  + 

Y  HEAD 

- 

"  P 

ESTAR 

e* 

E*_  in  Btu/Lbm 

K 

K 

Table  C-l 

KCM 

£cm 

Table  C-l 

KS 

Table  C-l 

KR 

K 

Table  C-l 

K-DISS 

Kd 

:  *0-Kcm-Kfi-Kr 

K  -  initial  kinetic 

£-Ln-fr.gy _ 
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■ROD  PROGRAM  MAINLINE 


C 

C 


_c 


.THIS  PROGRAM  COMPUTES  THE  BEHAVIOR  OF  A 
TWO-ELEMENT  ROD  KITH  A  DEFORMING  HEAD 

EXTERNAL, RODDE,OUTP 

COMMON  BMAX,5IGMA,ELOVO,ESTBK,LFREE,TH1CK,RHO,ESTRT,CD,FS 
_C OMMON  OBL1Q,  E S.T  A R ,  P P R 0 J , L B D 0 T , 1 L A V R ,_NLA_Y R ,L 8MF.T _ 


C 


.DIMENSION  PRMT (5) ,  Y ( 1 5) ,DERY ( l5)jAUKU/.1.5)_ 
DIMENSION  X(15), ITITL(40) 

DIMENSION  DEN3(5),ESTR(5),TH1K(5),DTAUI(5) 

DATA  NI  N,‘l  C  ASE»  NO  I  M/5 , 0  ,  15  / .  . 


fi 


1000  FORMAT  (1615) 

1001  FORMAT  (40A2) 

0j02 _ f  0RM4.T_<  st  1 0 1 0.) 


2001 
2002 
L.20Q3 
2000 
2005 


FORMAT 

FORMAT 

.FORMAT, 

FORMAT 

FORMAT 


(1H1 
UH 
(l  HO 
( 1 H  0 
( IH 


n 


,2M6_  FORMAT.. (1H0 
1  l 1 X, ' E* 

FORMAT  (1H 
FORMAT  OHO 
1  •  ESTR 1 3X 
FORMAT  oho 

1 


2007 

L200S 

2009 


P.  INC.') 

ROD  PROGRAM') 


54X, ' A.  R,  A. 

SIX, 'TWO  CELL 

40.A2.) _ _ _ 

JX,  'NCASE'  ,5X, 'NPRNT',5X, 'NLAYR'rRX, 'THICK') 
SX,I5,5X,I5,5X,I5,1X,E13,4) 

7.x , A R OD ». , 5 X ,  '.y.EL.OC U.yjL »9X,  »L/D' .  12X,  '  9 1 6_MA.’  »9Xf_'_E PS 01, 
FREE ', 1  OX, ' C  D ' , 7  X , • OBL 1 0 ' ) 

10X,6(E14.4,IX)) 

4X,  '..TARGET  LAYER' .,12, 6X, '.DENS',  1  OX  t  'THICK',  9X  , .  _ 

Till  •  ^ 


' DTAU ' ) 

7X,  'TAUMOX,  'VS*,  I  IX,  '  VCM« ,  10X,  'V« 


PERP',10X,'VT',12X, 


* 

2010  FORMAT  (IH 

1  7X, ' YHEAD' 

2011  FORMAT  (IH 

,  7X ,  •  PE N »‘,  1  0 X , MH  •  ^ V 1 X',' 
,6X, 'E8TAR*  ) 

jjWjJLK'»J2X,  'KCMl,lOX,il 

». 

2012 

2013 

2014 

FORMAT  (IH 
FORMAT  (IH 
FORMAT  ( 1H1 

, 9E14.4) 

) 

) 

2016 

FORMAT  (IH 

,1SX,7(E14.4,1X)) 

c* 

3001 

FORMAT  (IH 

, 'ERROR  4') 

L_3  0.0.2  _ 

F0RMA_L_UH5J  'ERROR  21L 

3004 

FORMAT  (IH 

, 'END  OF  LAYER', 12) 

c* 

3005 

FORMAT  (IH 

, 'END  OF  SHAFT') 

3006 

FORMAT  (IH 

, 'ERROR  3') 

c 

c 

c 

N0UT«6 

c 

READ*-; INPUT  CARDS' 

c* 

c 

10' 


-PISS') 


10' 


20 


Lil. 

c 

c 


.READ  ,  _ ___  _ 

ICASEa ICASEt 1 

IF  (ICASE,  LE.  NCASE)  GO  TO  20 

CALL  EXJT _ 

READ  (NIN,1001) 

READ  (N1N, 1002) 

READ.JNINt.lOOO) 

00  21  "n«T,nlayr 
READ  (NIN, 1002)  DENS (N) , ESTR (N) , THIK (N) ,DT AUI (N) 
CONTINUE _ _ _ ■_ 


ITITL 

VELOC, ELOVD,EFREE, SIGMA, EPSO, THICK, OBLIQ 
NPRNT,NLAYR _ _ _ 


INITIALIZE  X 


VECTOR 

1  17 


XU)  18  VS 
mi)  is  vcm 
X(3)  IS  VPERP 
X(R)  IS  VT 
X (5)  IS  MS 
X ( 6)  IS  MH 
X ( 7 }  IS  MS 
X ( 8 )  IS  B 
X  ( 9 ) _ IS  L 

x(io>  is  k 

XUD  IS  KS 

.  X ( 12)  18 _ KJCM _ 

X ( 1 3)  IS  KR 
X(IH)  IS  P 


C _ X  ( 15)  „1  S _V Hf ■AD _  _ _ _ _ _ _ _ _ _ 

c 

MXSTP«?000 

_ STOP*,  001  _  _ 

'  TQLER»,001  . . 

ZKSTRa.Ol 
ESTBKal . 0 

iKLB«r,0  . ~ . 

C0«1,0 

_  X ( 1 )«V6L0C  _  _ 

X  ( 4 ) '■  Z  *8  TR*V"E  IOC  ■  '  .  . . 

BZROaO.S/ELOVD 

ELZHO«ZKLB*BZRO 

I - E  LW8  a  EUZ  RB  7  b  2  R*0‘ - - - - - - - 

D)SS>ZKSTR*(4,0*EL0VR-C0,5m,*ELOVB*ELOVB/3'0)*ZKSTR) 

1  X(2)iVELOC-2,0*ELOVB*XC4) 

X (5)«0,0 

X(6)i8|0*ELOVD*ELOVD*BZRO*BZRO*ELZRO 

- x  ( 7 )  a  i  to  -k  (S') - - - 

X(8)«BZR0 
X (9) «ELZRO 

X~(  i'i)  aX ( 7) * X  (1)  *X’(  1 )  “ 

X(12)aX(6)*X(2)*X(2) 

XIaX(3).X(2) 

! - * ( 1  3 )  a x  (  6 )  *  <  175 * X  ( 4)  *  X  ( 4  )  ♦  X I  * X I ) 73 . 0 - 

X(10)aX(ll)tX(12)tX(l3) 

X(14)a0,0 

X( 15)a0»0  .  . 

DO  30  !■ 1 » NDIM 
VCI)aX(I) 

”30 - CONTINUE - - - 

C 

c 

t -  PR  I  NT  "THE "I  NPUT . 


C 

WRITE  (NOUT  #2001) 

- WRITE -(NOUT72O02)  .  . 

WRITE  (NOUT, 2003)  ITITl 
WRITE  (NOUT, 2004) 

- WRITE  ( NOUT i 2005)  NCASE,NPRNTiN7AVR, THICK 

WRITE  (NOUT, 2006) 

WRITE  (NOUT, 2007)  VE10C,ELCVD, SIGMA, EPS0,EFREE, CD, OBL1 

- DO~40~N177NL“A-YR  . 

WRITE  (NOUT, 2008)  N 

WRITE  (NOUT, 2016)  DENS (N ) , THIK (N) , E8TR (N) , DT AUI ( N) 

ire - CONTINUE 
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e  I 


,  <MWBWIlimiiiugi.ii  mih  mh. 


^KU'JtAinrc.t.K 


initialize  REMAINING  PARAMETERS 


„120_ 

c 


TAU«0.0 

TAUM»0,0 _ 

lLAYR«r 

LAYFGbO 

LCORRbO  _  _ 

LSHF  f  » 1 
FSaO.O 

09L1Q*OBLIQ*0,017«532R 
“'THENb*THfCK7C0S(0BlIG) 
OTAU«OT  AUI ( ILAYP) 

„RhO«.OeNJULAyP) _ _ 

ISTRTbESTRULAYR) 

TAULRbTHXK(IIAYR) 

B.MAXb  (1 1 0>£P_S ZRfl. _ 

2MOMO>VELQC 

zenobvelocaveloc 

KPT.-J _ _  _ 

ISTEPbI 

ITAUbO 

_ If  EXIT  BO _ _ _ 

LFLAGbO 

tTERMBO 

.„UPBNTbJ; _ 

UlNEBl243*(NLAYR«n 

LBDOTbI 

_ ®P_tPJ  55 _ 

DO  120  1  ■ 1 r  NDXM 
YU)bX(I) 

CONTINUE _ _ 


PREPARE  F ( 

DO  200  iSTEPi^iWlftF 

TAUbTAUM+DTAU 

PRMT(1)«TAUM 

PRMT (2) *T  AU 

PRMT(J)»DTAU/15, 

PRMT<<n«0,01 _ 

TEMPbI.O/NOIM 
DO  130  l* 1 # NDJM 

DERY(I)bTEMP _ 

CONTINUE 

INTEGRATE 


HPCG  INTEGRATION 


SYSTEM  OF  EQUATIONS 


C  _ 

"c  ~ 

c 


CALL  HPCG(PRMT»Y,OERY,NDIM,IHLF,ROODE»OUTP,AUX) 

IF  CIMLF.  IE.  10)  J?0__T0_lil _ _ 

WRITE  (NOUT #3002) 

CALL  EXIT 

fHt DK _F  Q R  T A Y  E R  'PEN E  T  R  A  T 1 0 N 

TAUOFBi|.0«TOL„ER),ATAULR  __  ___ 

IF  (YdSJt  LT.  TAUOK)  60  TO  135 
TAUBGbC1,OnTOLER)»TAULR 

IF  (YUi). J.E,_TAUBG)_60  T0_133 _ 

H  AT  10b  (  T  AUCR»X  ( 1’5 ) )  /  ( Y  (15 )  "XTfS )  ) 
DTAUbOTAU*RATIO 

TAUbTAUM _ _ _ 

GO  TO  100 
ILAYRbILAYRaI 
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OTAUiOUUlULAYR) 

RHO«DLNS( 1LAYR) 

ESTHT*ESTH(ILAYR) 

TAULP«tAULR»THlK(RAYR) 

7  '  LPRNTal 

1  L AYFG« 1 

IF  (1LAVR.  G7,  NLAYR)  LTERM* 1 
C 

C  STOP  THE  GROWTH  OF  THE  HEAD 

135  •  •  jjr“fVTiT)TTf*”vc i 4) )  V c iRjivus) 

IF  (Y(7),  6T.  0.0)  GO  TO  1359 

Y(7).«0*0 . . . „ . .  . 

Y ( 1 )a0.0 

Y ( 1 1 ) *0 .0 

1359  GO  TO  (1351, 1352, 136) , LBDOT 
*  1351  IF  Ct(8),  GT,  UMAX)  GO  TO  1352 
BPRED«2,0*Y(8)-X<8) 

IF  (BPREO,  LT.  BMAX)  GO  TO  136 
'  SLOPEa(BMAX-Y(6))/(BPRED-Y(8)) 

DTAU*(TAU«TAUH)* SLOPE 
LBD0T*2 

|  GO  TO  136 

1  1352  BMAXaY(B) 

OTAU«OTAUI (ILAYR) 

LBDOT  «3 
C 

C  PRE01CT  LAYER  PENETRATION 

rc  *  '  . . 

136  IF  (LAYFG,  EQ.  1)  GO  TO  1369 

YPRED«2,0*Y(15)-X(15) 

TAU0Kt(1.0-T0LER)*TAULR 
IF  ( YPRED,  LT.  TAUOK)  GO  TO  1369 
RATIOt(TAULR-Y(15))/(YPRED»Y<15)) 

f  . OT  AUP*R AT  10*  (T  ALI“T  .  ‘ . 

.  IF  (DTAUP,  GE.  DTAU)  GO  TO  1369 

DT  AU«DT  AUP 

IF  (LBDOT,  EQ.  2)  LBDOTil 
C 

C _ PREDICT^  ENp_ OF  SHAFT _ 

1369  GO  TO  (1370, 1371.1374). L9HFT 

1370  PDGMS«2.0*Y(7)-X(7) 

~  IF  (PDQM3,  GT .  '0,0)  60  TO  1373 

RATI0«Y(15)/(Y(I5)-PD0MS) 

DTAUP*RAT10*(TAU«TAUM) 

1 - 1 FCDTAUP."  fT~D T'A U )'  DT  A¥t DT  AUP 

LSHFT*2 
LPRNTal 

- GlTTOT3>3 - 

1371  L8MF  T«3 

•  Y(i)«0.0 

r - vnr>«r.u  ~  ~  “ 

v(ii)»o,o 

Y(4)«0,5*Y(8)*(Y(2)-Y(3))/Y(9) 

- rmOTiYH¥Ym 

IF  ( Y(4) .  GE,  FACT)  GO  TO  13715 
Y(4)»f  ACT 

- y  x  3  r«  y  ( 2POTm«v(  9)7Y(  s  r  ~~ 

13715  Y(13)«Y(6)*(0.5*Y(4)*Y(4)f(Y(3)-Y(2))**2/3.0) 
Y(10)aV(13)9Y(l2)9Y(l  1J _ 

TIME  STEP  AND  PREDICTIONS  FOR  NEXT  STEP  COMPLETE 
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z 

c 


c 

1373 


137  35 
13736 


1374 


13741 


r 

*  c 


J.J.7JL 


I  ISO 

c  " 

4  C 


L_c 

c 

<  C 


VB*Y(3)*2»0*(Y(2)»Y{3)) 

IF  (Vi,  GE,  Y ( 1 ) )  60  TO  13735 
F 6*1,0 

IF  (LCORR,  EO.  'I)  60  TO  13738 

Y(4)*0.5*Y(8)*((V(l>-VB)/Y(6)-(V(3)-V(2>)/Y(9)) 
_Y(13)«Y(6)*(Y(4)*Y(4)/2,0»(Y(3)*Y(2))**2/3.0) 
LCORRal 
GO  TO  13736 

F  S*0.0 _  _  __  .  _ 

LCORRsO 

CONTINUE 

..IF  U.8HF1,  GE.  3)  GO  TO  1374  ...  . 

IF  (Y(l),  6T,  VB,  OR,  Y(7),  GT,  0,5)  60  TO  1374 

lprnt*i 

GO  TO  1371  . .  _  .  . . 

TEMp»Y(5)-X($) 

IF  (TEMP,  EQ,  0,0)  60  TO  13741 
PPROJ«Y(14)U,950-Y(5)  )*  ( Y  ( 14 ) -X  ( 1 4 ) ) /T  EMp 
00  1375  Id'NDlM 
X  ( 1 )  a  Y  ( I ) 

..CONTINUE _ .  . . . .  .. 

TAUMaTAU 

_ _ CHECK.  FOR  RUN  TERMINATION  .... 

ZMOMaY (7 )* Y (1)+Y(6)*Y(2) 

RLAM0M«ZM0M/ZM0M0  „  _ _  __ 

“if  (ramom.  gt ,“ stop* “rso  to  ho 

LTERMal 

00  TO  ISO  _  ..... 

ZENER*Y(7)*Y(1)*Y(1)+YC6)*Y(2)*Y(2) 

raenh*zener/zeno 

IF  <RAENR.__GT,  STOP)  GO  TO  150  _  _ 

lVerm*i  '  . 

IF  ( Y ( |4) ,  LT,  THENO)  GO  TO  155 
LTCRM*! . . . . . . 

CHECK  THE  FLAGS 

“KPT.KPUl 

IF  (ISTEP,  EO,  l)  60  TO  190 
JF  (KPT.  NE.  NPRNT )  GO  TO  165  _  _ 

LPRNTal 

KPT*0 

TERM.  EQ.  1.  OR.  LEXI T , _E_0_*  1)  LPRNT*  1 _ 

IF  (LPRNT,  EQ,  0)  GO  TO  200~~ 

IF  (LFLAG.  EQ,  0)  GO  TO  195 


PRINT  THE  PAGE  HEADINGS 


IF  (LLINE,  NE,  0)  GO  TO  193 

UVl  TE  TfJ&UT  V20T4  ) - 

WRITE  (NOUT , 2003 )  ITITL 

WRITE  (NOUT ,2009) _ 

WRITE  (NOUT 1 2010) 

WRITE  (NOUT, 2011) 

_LL  IN  E  «LL  I  N£  ♦  4 _ 

L FLAG *6 


<  C 

c 


COMPUTE  AUXILIARY  QUANTITIES 


<  195  VB*Y(3H2.0*(Y(2)-Y(3)) 


wIDTH*2«0*Y (ft) 

TH1C*2.0»Y(9) 

KATM«1.0-Y(5) 

DISSK«VELOC*VEt.OOY(lll-Y(12)-Y(13) 

HE ADMa8,0*ELOVD*ELOVO* Y (8 ) * Y ( ft )* Y (9) 
RELK«Y(6)»(l,5*Y(«)*Y(<»)  +  (Y(3)-Y(2>)**2>/3.0 
EKRM*1.0-HEA0M/Y(6) 

IF  (ERRM.  IE.  .5,  OR.  LSHFT .  GE.  3)  GO  TO  1955 
WRITE  (NOUT . 3006) 

CALL  EXIT  __  _  _ 

F  1955  ERRK«l,'o«RELK/Y(13) 

E*C8TAR* 39. 91761 

i  C .  . .  ...  .  .  _  .  . . . 

C  PRINT  THE  OUTPUT 

C 

.  „  WRITE  (NOUT, 2012)  T AU, Y ( 1 ) , Y (2) , Y (3) , Y (4) , VB, WIDTH, THIC 

. WRITE  (NOUT,  2012)'  Y(14),Y  (6),  Y  (7  )‘,  Y  (5) ,  RATM,  Y  ( 1 5)  ,E  ' 

WRITE  (NOUT, 2012)  Y < 1 0) , Y ( 12) , Y ( 1 1 ) , Y ( 1 3) , DI8SK 
IF  (LAYFG,  EO,  0)  GO  TO  196 
ILAY«ILAYR»1 
WRITE  (NOUT, 3004)  IL AY 
_  LAYFG«0 _  __ 

196  ~  IF  (LSHFT ,  EO,  2)  WRITE  (NOUT ,30 05 ) 

WRITE  (NOUT, 2013) 

_ . _LLINE«LLINE*4 _ 

LPHNTaO  "  “  " . 

IF  (LLINE,  LE,  55)  GO  TO  197 
LFLAGH 

LLlNEaO  *  '  . . 

197  IF  (LTERM,  NE,  0)  GO  TO  10 

_ IF  (ISTEP,  EO,  1)  GO  TO  100  __ 

IF  (LEXIT,  EQ.  0)  GO  TO  200" 

IF  (ISTEP.  EO.  MXSTP)  WRITE  (NOUT, 3001) 

CALL  EXIT 

“*00  COKTTlNuE  . 

END 


RODDE— — ROD  PROGRAM  DERIVATIVE  EVALUATION 

SUBROUTINE  RODDEU# Y#DERY) 


THIS  SUBROUTINE  IN  THE  ROD  PROGRAM 
EVALUATES  DERIVATIVES  FOR  HPCG  INTEGRATION 

COMMON~BMAX, SIGMA, ELOVD,ESTRK,EF REE, THICK, RHO,ESTR 
COMMON  OBLIQ,ESTAR,PPROJ,LBOOT, ILAVR,NLAYR,LSHFT 

DIMENSION  DERY (15)# YUS) 


T,CD,FS 


lOfiO  ^PRMAI 
C 


LERROR- 


INITIALIZE  VARIABLES  FOR  THIS  STEP 


C..__ 


NOUT  *6 
VS«Y(1) 
VCM«.V<2)  ..... 
VPERP«Y(3) 
VT*Y(A) 

.EMb«i(5) _ 

EMH«Y(S) 

EMS*Y<7) 

_B*  Y_(5 ) _ 

EL*Y (9) 
EKiYUO) 

.EK3  F.YJUJ _ 

EKCMaVUZ) 
EKR«Y(13) 
PEN»YU<0  _ 
YHE AD»Y  US) 


PARAMETER  DEFINITIONS 


XUVPERP-VCM 

OELTA»VSiVPERPt2,0*XI 

»SO»B*B 

ELDSQbELOVD*ELOVD 


CHE  CK'  D13CRIMI  NAN  T  0  Fl  U  A  OR  AT  I C  T  E  R  M" 


BIGA««.O*B*EL/(2,0*8SQ-FS/i,0/ELOSa) 

BIGbiFS/(BSQ«ELDSQ*S,0-F8) 

OUADUC  3,042 (0*B1GA*RI6A)*EKR/EMH 

“00  A  62  ■IrCPflGS?  F  ATfifAC  f  ~ 

QUAD>0UAD1*QUA02 
IF  (OUAD,  GT,_0,O)__GO  TO  5 
WRITE  (NOUT, 1000) 

CALL  EXIT 

EVALUATE"  F  FACTORS 
IF  (LSHFJT,  LT,  3)  GO  TO  20 

ImsYoTo 

FS*0 • 0 

IF  (LSOOT,  LE.  2)  FH*OtO 
IF  CL  BOOT  ,  C(5,jrPH*l  ,0~ . 

EVALUATE  SEVERAL  DERIVATIVES 


IF  (F8«  EG,  0,0)  GO  TO  SO 


143 _ 


0C«YCn«»8ICMA/EM8 
GO  TO  60 
DERY<1)*0.0 

D£KY(5)«16,0*ELDSQ*B*EL*VT*FH 

DERY(6)*FS»DELTA-DERY(5) 

OERY(7)*-FS«DELTA 

DERY(U )sVS*(VS«OERY(?)42,0*EMS*ntRYtl)) 

0ERY(8)sVT*(1.0-FH) 

_ EVALUATE  EFFECT  I  VE_E*_  OF  J  hE_ROD  MATERI  AL 

IF  (EFREE .  EQ,  0.0)  GO  TO  65 
CALL  R00E8(VCM# VS#ESTRD) 

GO  TO  67 
ESTRD«0.0 

_  ~-At  j  BACkF ACE  EFFECT  .  *  . 

EST  AR*ESTRT 

IF  (HAY R V  L f "«~~NL A  YR )~GO  TO  90 

BRKPT«THICK-R.O*B 

PENTH*YHEAD*COS(OflLIQ) 

'IF' (PENT H.  IE.  THICK)  GO  'TO "70  . ~ 

ESTARaO .0 
GO  TO  90 

IF  (PENTH.'VE.  B'RKP'tTGO  to  90 . . 

E3TAR«ESTRT*(),0-(PENTH-BRKPT)/R.0/B) 

IF  (YHEAO,  LT.  PEN)  GO  TO  93 

IF"  ( VPERP .* LT 0 .  O')  GO  TO  93  .  . 

force*estar>co*v»erp*vperp/2.o 

OER Y ( 1 U )» VPERP 

GO'TO  95 . 

FORCEiO.O 
DER Y ( 1 9 ) ■ 0  •  0 

ETAiVT - - -  -  " 

X1*VPERP-VCM 

ETA2«ETA*£TA 

XI2«XI*XI  . . . . 

ZKRl»«DERY(5)*FH*(ETA2+XI2/3.0) 

ZKR2*«8.0»RHO*BSQ*EL08Q*XI*FORCE 

fact«vs-vcm  . . . 

ZKR3«2.0*SIGMA#F8*FACT 

ZKRU*»2.0*l)ERY(7)*(FACT*FACT/2,0"E3TRD) 
f)ERY(  U)«ZKRI4ZKR24ZKR34ZKRR 

DERY(2)*-VCM4(DERY(5)40ERY(6))-VS*DERY(7)-EMS«0ERY(I) 

I  -<l,O*BS0*ELO3Q*FORCE*RHO 

0ERY(2)*DERY(2)VEMH  . *  ‘ 

DERY(12)«VCM*(VCM*OERY(6)42,0*EMH«DERYC2)) 

OERY(10)»DERY(il)4DERY(|2)4DERYU3) 

DERY(9)»VPF.RP-VCM 
OERY ( 15)»VPERP 

.  FACTORS  FOR  VPt  HP  AND  VT“ DERI  VAT  IVES  . ' 

ABAR*(3,n*0ERY(13)-(1.5*ETA24Xl2)*DERYC6))/EMH 
Cii(OERVti)-0ERY(2))tF3/EMH 
FACU«DERY<a)*VT/B 
F ACT2«XX*(DERY<8)4VT)/EL 

C  B  A  R  ■  C'  l-2,'0*(  FAC  T  1+  FA  C  T  2  )7B  ‘  . 

FACTUDERY(6)*(2.0*VT/B4XI/EL)/EMH 

FACT2iXl2/EL/EL42.0*VT*DERY(8)/B3U 

CBAR*CT-FACTf4FACT2 

DB ARb 1 . 0/B  144 


EBAHU.0/EL-F3/EMM 


LAST  TWO  DERIVATIVES 
DEN0Ma«70*Xi*DBAR-5,6nETA*EBAR 

DERY(S)#DERY(2)Y  (2, ABAR*DB AR-3, 0*ET  A*CBAR) /DENOM 

DERYi4)s(CBAR.EBARo(DERV(S)«DERY(2)))/a.0/DeAR 

RETURN 


uuuu 


c  Roots— —ROD  PROGRAM  ESTAR  OF  ROD 
SUBROUTINE  R0DES< VC  1 , VS1 ,ESTR) 

THIS  SUBROUTINE  IN  THE  ROD  PROGRAM 
CO-VUTES  THE  EFFECTIVE  ESTAR  OF  THE  ROD 

COMMON“BMAX, SIGMA,  ELnVD,ESTbK,tFREE,T HI CK,RHO,ESTRT»CD»F3 
COMMON  OBLIQ,ESTAR,PPROJ,LBOOT, I L A YR , NL AYR , l UHF T 

R  AT  I  OVoTsTcTsl^Vrrr*  *  Si/E  FREE 
IF  (RATIO.  GE.  ESTBK)  GO  TO  10 
tST»«EFREE*RAT10/fcSTRK  _ 

RETURN 
10  ESTR*EFREE 

_  RETURN  _ _ _ 

'  END 


TTs*c"rarn0B3ES! 


IS  *-r^;  ,. 


r-JVKEHsvi'3  T.zxBssrraFnnxnf’? 

. .  «  •  -*-  *  «*wm  •ir.'r;W3*J«eW3  rrV-  /¥*.***>*, 


SUBR0U1 INE  OUTM(Xr YrDERT#IHt.F,NOIHfPRMT) 
RETURN 
r  wn 


■■■{yj  r-ia  i-.Jw'.agariJSfit'j  a-* 


T^.nwr-wfrm»rMgw  'jH^^gg^w,  vxam-  ■ 


MFC  6  -  PREDICTOR-CORRECTOR  ORDI NARY  DIFFERENTIAL  EQUATION  SOLVER 
SUBROUTINE  HPCC (PRMT , Y , DERY , NDI 1 # IHLI # FCT # OUTP# AUX) 

DIMENSION  PRMT(5)tYU)rDERV(l)fAUX(!6rl) 

NDIMiNDIl 
Ns  ] 

lHLFsO 
XsPRMT ( 1 ) 

HsPRMT ( 3 ) 

prmt;5)«o,o  __  _  _  _  _ _ _ _ 

DO  T  i»i#NOTM  ■"  " . .  *  '  ~ . 

AUX<16#I>*0,0 

AUX(15#I)*DERY<I> 

AUX(t,I)*Y(l> 

IF  IH*(PRMT(2)-X))  5,2, A 

IHLFS12 

60  TO  4 

IHLF-15 

CALL  FCT(X,Y,OERY>  „  . 

CALL  OUTP (X,Y, DERY, IHLF,NDIM,PRMT) 

IF  tPRMT(S))  6,5,6 

_ IF  (IHLF)  7,7,6  _ _ _ _ _ 

IHLIsiHLF  . . 

RETURN 

_D0  8  I*1,NDIM  __  __  _ 

AUX(8i,I)«DERY(i  j 
ISWsl 

_60  TO  100 _ _  _  _ 

XsXVh  "  ■’ 

DO  10  I*t,NDIM 
l  AUX  ( 2 , 1 J)  *Y  C 1 ) 

IMLFslHLFVl  ’ 

XcX-H 

_  DO  12  I *1 , NDI M _ _ _ _ _ 

f  “luriTTTTi'Aux  ("271 )' 

H»0,5*H 

N*1 


lSWsj 
60  TO  100 
X*X*H 

- C~A  L  CT f T  TTf7Y7  DE  R  Y ) - 

Ni  2 

DO  14  Ul'NDlM 


AUX (£ , 1 ) s Y (  I ) 

AUX (9 , 1 ) sDERY  ( I ) 

2  SN*S 

- 6o~tO~TC6 - 

DELT«0.0 
DO  16  1*1 »NDIM 

DELT *DELT+AUX <15,11  *ABS(Y ( I J -AUX  (4,1)1 

DELT*0,06666667*DELT 

IF  (DEL7-PRMT (4) )  19,19,17 

— rrrfHLT-roTTi/m  re - 

IHLF*11 

XsXfH 


X«X  <H 

CALL  FCT(X, Y»  DERY) 
— DO  '  2(T  1THTID 1 « 
AUX(3«1)*YCI) 

«UX ( 1 0, I) *DERY( I ) 


2.3 _ _ 


*  soo 
i  '  !~2o  r 

-■)  202 


GO  TO  100 

Nil 

XaX4H 

CALL  FCl(X,YrDERY) 

XtPRMT ( 1 ) 

00  22  1*1, NDI* 

AUX(11,I)*DERV<1) 

Y(I)*AUX(l, 1 )4H*< 0.375* AUX (8,1) ♦0,791 6667* AUX<9, 1) 

1  -0,2083333* AUX (10, I ) 40.04 166667*DERY ( 1 ) ) 

_XsX4H _ _  ..._  _  _ _ _ _ 

N«N4  i 

CALL  FCT(X,Y,OERY) 

_C  ALL  QU1P(X,Y,PERY,1HLF,NDIM,PRMT) _ 

IF  (PRMT ( 5) )  6,24,6 
IF  (N-4)  25,200,200 

„0.0.26.I«i,.NDIH _ _ _ _ _ _ _ _ _ 

AUXIN, X ) *Y ( I ) 

AUX (N47, I )*DERY(I) 

„1F(N-S)  27,29,200 _ _ _ _ _  _ 

DO  28  I*1,NDIM 
DELT*AUX(9,I)4AUX(9,I) 

_DELj«OEL!tDEL.T _ _ _ 

Y(I)*AUX(1,I)40,3333333*H*(AUX(8»I)4DELT4AUX(10,X)) 

GO  TO  23 

_.DO_30_  I  *  1  ,.ND  1 M _ 

DELTaAUX (9, X)4AUX(|0,1) 

DELT*DFLT4DELT4DELT 

.  JCJ I }  >AUX.(JjlS  )-±-C .  JXS  *H*J  A  U  XHb.  UJt  &l  LI  ±A  U  X_(  .1  X/J.U _ 

GO  TO  23 
DO  101  I»1 , NDJM 

.2*H*AUX(N47jD- - _  _ _ _ _ _ _ _  .... 

AUX(5,I)«2 
Y(I)*AUX(N, 1)40, 4*2 

_2*XiP,4*H _ _ _  _ _ 

CALL  FCT(Z, Y,DERY) 

DO  102  1*1,ND1M 

_z«H*DER.yj.n _ _ _ _ _ 

AUX(6,X)*Z 

Y(1)*AUX(N, I )40,2969776*AUX(5, 11*0,1587596*2 

Z«X40.45S7372*H _ 

CALL  FCT(Z,Y,DERY) 

DO  103  I«l,NDIM 

__2«H*DEmiJ _ _ - 

AUX(7,n*Z 

Y(1)bAUX(N, 1)40, 2181 004*AUX(5,X )-3,050965* AUX(6,1)43,832865*Z 

_Z»X4M _ 

CALL  FCT(2,Y,DERV) 

DO  104  I«i,NDIM 

„JlU)aAUX<N,I)40,1747603*AUX(5,I)-0,5Sl«B07*AUX(6,J>  _ 

1  4l,205536*AUX(7,l)40o17118«B*M*DERY(l) 

GO  TO  (9, 13,15,21), ISW 

_  X3TEP*3  _ _ _ _ _ _ 

IF* (N -8)  204,202,204 
DO  203  N*2, 7 

_ DO  203  1«1,ND1M _ _ _ _ 

AUX (N-i ,I)*AUX(N,1) 

AUX(N46,1)«AUX(N47,I) 

_N*7 _ 

N*N4  1 

DO  205  I«1,NP1M 

AUX  (N-l ,  I  )»Y  ( 1 )_ _ _ _ _ _ 

AUX (N 46, i)*OERY  (1 )  ,.Q 

X*X4H 


201 


213 

210 

_ 2.1 5 

216 

217 

p21« 


_ 


6  ISTbPslSTEPtH 
00  207  1*1 iNOIM 

*  •JSKJSJiriJ  siiJisJJ????*"*  «•**<•«*'  * »—» o'**. '  i  > 

V(l)*DtLT-0,9256l98*AUX(16,T)  '  - 

7  AUXCl6,l)«DtLT 

...  CALL  FCT(X,Y,DERY) 

00  206  1*1, NDIM  - 

1  *  •*  .-H.CMII** .  > >  t«**. , , 

AUX(16,I j*AUX(16,I)-DELT 
i  Y(I)sDELT+0.0?«38017*AUX(16,I) 

...  .  DELTstO  ,  0 

00  209  1*1, NDIM . 

>  0ELT*DELT*AUX(15,1)*ABS(AUX(16,1)) 

_ IF.fpELT-PRMTlO))  210,222,222 

>  CALL  FCT  (X#  Y,DER~Y) - - - - - 

CALL  OUT P(X,Y,0£RY,IHLF,  NDIM, PRMT) 

_ 1  P_  (P W.LC5U.  21 2# 2 11 , 21 2 

IF  (IHLF-11)  213,212,212  - - 

!  IHLI*IHLF 

RETURN  _  _  __ 

'  I F  tWt X-PRMT ( 2) ) )"2iIV 21*2 , 2T2 
i  IF  (ARS(X*PRMT(2))-0.1*ABS(H))  212,215,215 

_ IF_<DELT-J.02*PRMT  («))  216,216,201 

IF  ( IHLF )  201,201,217  ‘  ~  - 

IF  (N-7)  201,216,218 
_ IfL.USTEP-t}  201,219,219 

I  MOO* I S  T  EP/ 2  - - - 

IF  (ISTEP-IMOO-IMOO)  201,220,201 
H*H4H 

IHLF*IHLF«i  - - - - - - - — 

1STEP*0 

DO  221  1*1 , NDIM 

TuX  i  - - 

AUX(N-2,I)*AUX{N-9,I) 

_ AUX(N-3£I)*AUX(N-6,I) 

AUX (Ny6, X )*AUX(N*5,  X)  . —  -  -  - . — 

AUXCN45,I)*AUX(NO,I) 

_ AUX(NtO,I)*AUX(N  +  l,I) 

CELT* AUX (Nt 6 , XT' ♦  AU  X ( N*5, I ) 

OELT*DELT*DELT+DELT 

lA4AuX(NjJ*ll’-2963t(V  -  ,"AUX<N"3'I,,*3,3611J^~  tCEP7CIl4DELL_ 

60  TO  201  * 

IHLFbIHLFYI 

IF  ~riHCF^lDT~2«T223r2  - - 

H*0,5«H 
_ 1 8  T  E  P  *0 

00  229  T*1 ,  NDIM  . - - -  -  -  - 

1  v!i.,)M?w023?!?25*f?S:2«u?<N"1,1>*l35*0*AUX<N,2'I,+floto*Aux(N.3,i) 

~2~^AWt  — 1  ® *  — ( Nt6>I^fe.0*AUX(N45,I) _ 

/!JfJSI*!l:i.,,!?:023??625*(12*0*AliX(N*,»I)4j:»5.oAAux(N-2,n 

-i— I> > -0. 0239375* (AUX(Nt6, I) 

2  ♦r8,0iAUX(NY5,I)-9.0AAUX(NY971J)«H— - 

AUX(N-3,I)*AUX(N.2,n 

AUX(NY9,I)*AUXCNfS,I) 

~ X*X-H - - - - - 

OELT«X-(HYH) 

CALL  FCT(DELT,Y,OERV) 

— 0(T2^5~r*  17N  DIM - — - 

AUXIN-2, I)*Y(I>  150 


225 


AUX(N+5,I)sDERY(Z) 

Y<I)tAUX<N-<l,!) 

DELT«D£I.T-(NtH) 


CALL_FCHDELT#t,DERY) 

DO  226  l«Y*NOXM  "  '  " 

DELT*AUXCN*5#IM  All  X(N  ♦<!»!) 

0ELT«DELl40EL7^DEUT  _ 

AUX(16,  n*e.962963>(AUX(N-l,  I)-Y(I)  j-3. 3*1 U 1«H« ( AUX (N+6r I ) 


} 

rt 


1  4DELTfOERY( I ) ) 

226  AUX(N+3,l)«DERY(i) 
60  TO  206 
END/ 


